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Zodlogists have long known that the purine body guanine, alone or in 
some combination, is responsible for the silvery appearance of areas of skin 
in most fishes. It is also a familiar fact that crystals of this substance are 
contained in special cells, variously called guanophores, leucophores or 
iridocytes, which respond, as do other chromatophores, to stimuli received 
through the eyes. The paling of a fish upon a white background results 
in part from the withdrawal from view of dispersed melanin granules 
within the melanophores, but also in part from the centrifugal movement 
of guanine particles within the guanophores. Reverse relations, with re- 
spect to both substances, are manifest upon a black background. 

That the same stimuli which call forth these rapid adjustments may 
result, when long continued, in actual large-scale changes in the amount of 
pigment present has been shown by a number of workers.' For melanin, 
fairly exact relations have now been established in certain fishes between 
the amount of this pigment and the albedo of the background to which 
they have been subjected. 

For guanine, no exact relations of this sort have hitherto been deter- 
mined, though several workers have recorded the occurrence of an increase 
in the number of guanophores or in the quantity of visible guanine in ani- 
mals kept in white containers.? In general, however, guanophores have re- 
ceived relatively little attention from students of chromatophore activity. 
This is doubtless largely due to their inconspicuousness in comparison 
with melanophores, likewise to their ordinarily being rendered invisible 
in histological preparations. The only worker known to me who professes 
to have undertaken the measurement of guanine in connection with color 
changes in animals (Meyer, op. cit.) reports a purely negative result, con- 
cluding that any possible experimental differences were swamped by the 
individual variability among the animals studied. This writer fails to state 
the method which she employed in making these measurements. 
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After the somewhat extended assays which have been made in this labora- 
tory of the black pigment melanin, in relation to the visual environment, 
a comparable series of quantitative determinations of the white pigment 
guanine seemed called for.* 

A local marine fish, Girella nigricans, was employed for this purpose. 
Two hundred and forty sexually immature* specimens were apportioned 
among 20 painted glass bowls.’ Of the latter, there were four each of five 
shades, ranging as follows in respect to albedo, relative to white as a stand- 
ard. 
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IS SS é ya awas ons Vary § 4nd dale 9 RS pas 17.80 
TIE OO ELEC ES PRES 6.90 
eign a dino cath bakkie ws bus Kaen eures 1.42 


Each set of four nearly identical bowls was kept in a separate cabinet, 
lighted by a 100-watt lamp overhead. Running sea-water was supplied 
to the bowls. 

Of the original number of fishes, 217 were alive at the end of four months, 
when they were killed with chloroform vapor, dried with paper towels and 
weighed, commonly in lots of six or seven. These were then dipped into 
boiling distilled water for three-fourths of a minute and the “skins,” in- 
cluding the scales and fins, were removed and dried for several hours in an 
oven at 100°C. These dried ‘‘skins’’ were also weighed to the nearest 
tenth gram. 

The fishes from the first bowl in each cabinet, which were fewer in num- 
ber and much larger than the others, were sacrificed for preliminary tests of 
method. The occupants of the other three bowls of each set were divided 
into two lots of six or seven fishes each. There were thus six lots from each 
cabinet (thirty lots in all) which were available for guanine assay. Of these 
six lots only five have been used. 

Thus the fishes under consideration are designated as follows: ‘‘1,” 
2.1,” “2.2,” “3.1,” “3.2,” “4.1,” “4.2”. Of these the ‘1’ lots were used 
for preliminary tests, while the ‘‘4.2’’ lots have not yet been used at all. 

The matter is further complicated by the fact that the solid material of 
the original ‘‘2.1” series was divided into two equal portions (‘‘2.la” and 
‘*2.1b”’) and each of these was treated as an independent sample. The re- 
sults from these are not quite comparable with those from the subsequent 
series, in that the extraction with acid from the former was less complete 
than from the latter. 

For the methods which I have employed in the extraction and measure- 
ment of guanine I am indebted to Dr. George H. Hitchings of the Wellcome 
Research Laboratories, Tuckahoe, New York, who is reporting upon his 
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own studies of some of the same material in another article in this issue of 
the PROCEEDINGS.’ It is a pleasure to thank Dr. Hitchings for help and 
advice, kindly given at various stages of this undertaking, though I cannot 
hold him responsible for every minor step in my procedure. I must also 
acknowledge the important help of my colleague, Dr. Denis L. Fox, whom 
I have consulted frequently throughout these studies, as well as using some 
of the facilities of his laboratory. I must mention, too, the valuable tech- 
nical assistance of Miss Ruth A. Larson during the earlier stages of these 
investigations. 
The procedure here adopted was as follows: 


1. The oven-dried skin was extracted with V/1 H2SOx,, in the proportion 
of 10 cc. of the dilute acid per gram of the original wet weight of the fishes. 
Extraction here comprised grinding of the dampened material in a mortar, 
followed by four cycles’ of prolonged shaking and centrifuging, in alter- 
nation, using a new fraction of the acid for each cycle. 

2. Measured portions of this acid extract were heated for one-half hour 
in 100-cc. centrifuge tubes in a water bath. 

3. These were made slightly alkaline with NaOH, using phenolphtha- 
lein as an indicator. 

4. The contained purines were precipitated with copper sulfate, fol- 
lowed by sodium bisulfite, at boiling temperature. 

5. The precipitate was brought down by centrifuge. 

6. Precipitate washed twice with boiling water, followed each time by 
centrifuging. 

7. Precipitate redissolved in boiling dilute HCl. 

8. Copper precipitated by H:2S, passed through the heated solution. 

9. Last precipitate removed from the guanine-HCl solution by filtra- 
tion through a sintered glass filter. 

10. Phenol reagent* (Folin-Denis-Ciocalteu) added to the filtrate, 
followed by NazCOs3. 

11. Colorimetric determinations made of the resulting blue solution 
by means of a Bausch and Lomb Visual Spectrophotometer. 


The degree of specificity to be expected from this method is discussed by 
Hitchings and Falco (pp. 294-297 of this issue). These workers have likewise 
added independent evidence, by the application of an altogether different 
method, that we are actually dealing here with guanine. This is impor- 
tant, since the procedure which I have employed is not sufficiently specific 
to distinguish between guanine and certain closely related compounds. It 
is gratifying to note that the values obtained by Hitchings and Falco, by 
the use of these two methods agree very closely; and it is further gratifying 
to note the reasonably close agreement between their figures and mine for 
the four samples which were.assayed independently in the two laboratories 
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by the copper precipitation method (Black 3.1, 3.2; White 3.1, 3.2). My 
figures should be compared with those in the middle column of Hitchings 
and Falco (0.285, etc.). 


In making the spectrophotometer determinations, absorption values 
(logarithmic) were recorded by me at 550, 600 and 650 my. Colorimetric 
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FIGURE 1 





Readings with dilutions of standard guanine-HCl solutions. Abscissae = spectro- 
photometer readings (logarithms of absorption at 550, 600 and 650 my). Ordinates = 
mg. guanine per 100 cc. of solvent. 


analysis of this sort obviously requires the previous plotting of graphs, 
based upon various dilutions of a standard guanine solution (Fig. 1). The 
values given in the tables, and in figures 2 and 3, are the average values ob- 
tained for these three wave-lengths. Five readings at each wave-length 
were comprised in each determination. It must be added that two to five 
subsamples of the acid extract were employed in the case of every lot of 
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fishes, and that the means here presented are based upon all the values thus 
obtained, commencing with ‘‘2.la,’’ except for a few cases in which the 
treatment had been varied for special purposes. 
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FIGURE 2 
Guanine values for Girella series. Abscissae = albedos of the five shades of 


background. Ordinates = mg. guanine per g. of wet weight of fishes. Continuous 
line based on means of series ‘‘2.2,’’ ‘‘3.1,” “3.2,” ‘‘4.1.”" Range of each series indi- 
cated by vertical bars. Broken line based upon means of ‘‘2.1a”’ and ‘‘2.1b.” 


TABLE 1 
‘GUANINE, First EXTRACTION 


BLACK DARK GRAY MEDIUM GRAY PALE GRAY WHITE 
mc./G mG./G MG./G MG./G mc./c MG./G MG./G MGc./G MG./G MG./G 
TOTAL DRY TOTAL DRY TOTAL DRY TOTAL DRY TOTAL DRY 


wT. SKIN wT. SKIN wT. SKIN wT. SKIN wT. SKIN 
2.1a 0.175 0.243 0.366 0.604 1.036 
2.1b 0.179 0.264 0.397 0.618 1.021 
Means 0.177 0.253 0.381 0.611 1.028 
2.2 0.231 3.61 0.286 4.27 0.3871 5.80 0.756 10.65 0.972 13.89 
3.1 0.260 3.71 0.293 4.58 0.483 6.90 0.766 11.27 1.084 15.27 
3.2 0.239 3.62 0.300 4.48 0.504 6.00 0.851 11.35 0.913 13.68 
4.1 0.256 3.77 0.333 5.05 0.467 6.97 0.851 12.33 1.110 16.32 
Means 0.246 3.68 0.303 4.59 0.456 6.42 0.806 11.40 1.020 14.78 


While occasional considerable differences were encountered between 
subsamples which should have yielded identical results, the mean deviation 
of such values from their average was 2.5 per cent (extremes: 0 and 6.6 
per cent). Such differences are mainly small in comparison with the differ- 
ences among the various experimental lots which had been subjected to dif- 
ferent backgrounds. 
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The mean guanine values of the various lots of fishes herein discussed are 
presented in table 1. These values have been computed in two ways. In 
one column.-is presented the amount of guanine in the ‘‘skin’’ (as defined 
above) per gram of total original wet weight of the fishes. This may not seem 
to be a very logical value to employ, since the very considerable amount 
of guanine present in some other tissues of the body was not ordinarily 
known. In the other column we have the amount of guanine in the skin, 
in proportion to the dry weight of the skin. This last would doubtless be the 
preferable value to employ, except for the circumstance that the ‘‘skin,”’ 
as here dealt with, was inevitably accompanied by a small but rather vari- 
able amount of adhering muscle. This would somewhat affect its 
weight, probably without correspondingly affecting its absolute guanine 
content. Thus the values in the first of these columns seem to furnish a 
fairer index of the changes in the concentration of this substance. Actually, 
the figures based upon the dry weight of the skin yield a graph very similar 
to those shown in figure 2. 

No calculation of probabilities seems necessary in support of the rela- 
tions shown by these graphs. Considering the four fully comparable 
series, ‘2.2,”’ “3.1,” “3.2,” and “4.1,” we find that the guanine content 
of fishes from the white bowls is 4.15 times that of those from the black 
ones, when ratios to total wet weight are considered; 4.02 when ratios to 
dry skin weight are considered. 

In view of the readily visible guanine deposits in the peritoneum of these 
fishes, it is not surprising that considerable of this substance may be ex- 
tracted from their “bodies” after removal of the “skins.” Only two lots 
were tested for this purpose, but in these the total amount contained in the 
“bodies” was absolutely greater than that in the ‘‘skins” (1.06 and 1.19:1). 
Since the latter constituted only about 23 per cent’ of the dry weight of the 
fishes, the concentration of guanine was obviously much greater in the 
skins. 

It was early found that even after four ‘‘cycles’’ of extraction, as above 
described, the fish material continued to yield some substance which gave 
distinct traces of the blue reaction with phenol reagent and sodium carbon- 
ate. The material from the ‘‘4.1”’ series was subjected to these additional 
extractions, these being carried through four cycles, as had been done 
originally. The results are shown in table 2. 

The mean absolute values of these ‘‘second extractions” thus do not differ 
significantly for the different lots. Indeed, their magnitude appears to 
bear no relation to the source of the chromogen, a circumstance which is in 
striking contrast to the figures for the “first extractions,’ which range from 
0.256 to 1.110, as we pass from black to white. This fact suggests that we 
are dealing with a different chromogen in the earlier and later extractions, 
and that in these last the amount is not influenced by optical conditions.” 








as 
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I may add that a third extraction (i.e., series of extractions) yielded 0.009 
and 0.008, respectively, for ‘‘black’’ and ‘‘white’’ lots (the other lots were 
not tested). These last figures are of about the same magnitude as was 
obtained when the H2SO, alone, of the same dilution, was employed, with- 
out any known contamination by organic matter. It would seem that all 
of these guanine measurements are subject to a small constant error, which 
perhaps could be determined more precisely and applied to the results pre- 
sented in the tables. 

Several earlier writers have stressed the reciprocating relations between 
guanophores and melanophores, resulting from changes in the background, 
but without attempting to measure either of the pigments concerned. Thus 
Kuntz (1917) remarks that ‘complete adaptation to a white background 
involves a rearrangement and probably an increase in the number of the 
guanophores and a material reduction in the number of melanophores in 
the layer just beneath the epidermis” (op. cit.). 


TABLE 2 


GUANINE, SECOND EXTRACTIONS 


FIRST SECOND EXTRACTIONS % OF FIRST 

EXTRACTIONS 1 2 MEAN EXTRACTION 
Black 0.256 0.021 0.017 0.019 7.4 
Dark gray 0°. 333 0.024 0.020 0.022 6.6 
Medium gray 0.467 0.024 0.024 0.024 5.1 
Pale gray 0.851 0.022 0.019 0.021 2.5 
White £110 0.025 0.018 0.021 1.9 


Murisier (1920-1921) notes that melanin and guanine seem to be sub- 
stances which may substitute for one another, the amount of one which is 
formed being inverse to that of the other. Like melanin, it is formed under 
the influence of retinal stimuli. Murisier suggests that both guanophores 
and melanophores form an accessory excretory system, but one which re- 
tains its products instead of discharging them. 

Millot (1923) sets forth these same ideas at greater length: “The three 
pigments (melanin, guanine, lipochrome) have almost the same physio- 
logical value and belong to the category of waste materials’ (trans.). Al- 
though guanine is an excretory product, its abundance in the tissues is al- 
together independent, he asserts, of the food consumed or of the activity 
or integrity of the renal apparatus. 

“Experiment reveals in a number of cases the existence of an equilibrium 
between the melanin and the guanine such that any diminution of the 
quantity of one pigment brings about a compensatory augmentation of the 
other.”’ 

While this reciprocal relation between guanine and melanin has been 
affirmed by several previous writers, it has never, so far as I know, been sup- 
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ported by quantitative data. Although Millot discusses a method for iso- 
lating guanine, he gives no figures to support his claims that the amount of 
this substance may be influenced by visual stimuli. 

In figure 3 of the present paper I have placed in comparison the curve for 
guanine production in Girella nigricans with that for melanin production 
in this same species, which was obtained under almost identical experimen- 
tal conditions.‘' The approximate character of the values upon which 
both of these curves are based has been freely admitted. 
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Comparison of curves for guanine and melanin in Girella. Guanine curve = con- 
tinuous line of figure 2. Melanin curve is taken from Sumner, 1948. 


While the melanin values, from black to pale gray, vary inversely as the 
- logarithm of the albedo, no such arrangement holds for the guanine values. 
However, the two curves agree in displaying a marked change in their 
trend beyond pale gray, that for melanin ceasing to fall (indeed, rising 
slightly) at this point, that for guanine undergoing a marked reduction in 
the rate of increase. It is, of course, unfortunate that an intermediate al- 
bedo between “‘pale gray” and white was not intercalated in this series, but 
the relative length of this interval was not realized by eye at the time of 
mixing the paints. 
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The outstanding fact revealed by these experiments is that the course of pro- 
tein metabolism, and the production of two important metabolites, melanin and 
guanine, are strongly influenced by the visual field of the fish, and influenced 
in such a way as to make for the animal's concealment. 


These fundamental organic processes are thus caused to change their course 
in the interest of superficial appearance, and this with no apparent effect upon 
the growth of the fishes. It is true that the strict accuracy of this last state- 
ment should be tested more carefully than was done in the present experi- 
ments, in which the original mode of selection of the animals, at the time 
they were distributed to the various bowls, appears to have led to minor 
size differences which are irrelevant to the question at hand. 


That this mechanism for concealment may be of life-saving value has been 
shown by considerable experimental evidence.’ 


If these two pigments (three, if we follow Maillot) actually function as excre- 
tory products, it is surely surprising that the necessary excretory balance should 
be maintained through reciprocal adjustments in the production of substances 
so widely unrelated chemically. Plainly, there are physiological (shall we not 
say psychobiological?) problems here which call for extensive future investiga- 
tion. 


* Contributions from the Scripps Institution of Oceanography, New Series, no. 235, 

1A summary of observations in this field was published by Sumner, Biol. Rev., 15, 
351-375 (1940). Some later results have been reported by Dawes, Jour. Exper. Biol. 
18, 26-49 (1941); by Sumner and Doudoroff, Biol. Bull., 84, 187-194 (1943); and by 
Sumner, Biol. Bull., 84, 195-205 (1943). 

2 Kuntz, Bull. U. S. Bur. Fisheries, 35, 1-29 (1915-1916); Murisier, Revue Suisse de 
Zoologie, 28, 45-97, 149-195, 243-299 (1920-1921); Millot, Bull. Biol. de la France et de 
la Belge, 57, 261-363 (1923); Hewer, Phil. Trans. Roy. Soc. London, B215, 177-200 
(1927); Meyer, Zool. Jahrb. (Abt. f. allgem. Zool.), 49, 231-270 (1931); Sumner and 
Wells, Jour. Exper. Zool., 64, 377-403 (1933); Odiorne, Proc. Nat. Acad. Sci., 19, 750—- 
754 (1933). 

3 | think that guanine may here be legitimately rated as a pigment, inasmuch as it 
actually functions as a pigment, i.e., influences the color (in the broad sense) of the ani- 
mal. 

4 This point is important to note, since it makes certain that no sexual differences in 
pigmentation can have complicated the picture. 

5 This is the same set of bowls that was employed in the writer’s melanin experiments 
with the same species (Biol. Bull., 84, 197 (1943)). They were of clear, colorless glass, 
painted on the outside, and sand-blasted within to eliminate gleams of reflected light. 

6 Hitchings and Falco, Proc. Nat. Acad. Sci., 30, 294-297 (1944). See also Hitchings, 
Jour. Biol. Chem., 139, 843-854 (1941); Hitchings and Fiske, 140, 491-499 (1941). 

7 Only one cycle in the case of ‘‘2.1@” and ‘‘2.1).” 

8 Folin, Laboratory Manual of Biochemistry, D. Appleton-Century, 1934; Folin and 
Ciocalteu, Jour. Biol. Chem., '73, 629 (1927). 

® Based upon only one lot of 3 fishes. 

10 Of course, a certain amount of this second chromogen may be supposed to occur 
in the earlier extractions. 
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11 The latter from Sumner, Biol. Bull., 84, 195-205 (1948). 
12 Sumner, Amer. Nat., 69, 245-266 (1935); Proc. Nat. Acad. Sci., 20, 559-564 (1934); 
21, 345-353 (1935). 


THE IDENTIFICATION OF GUANINE IN EXTRACTS OF 
GIRELLA NIGRICANS. THE SPECIFICITY OF GUANASE* 


By GEorGE H. HITCHINGS AND ELvira A. FALCOo 
THE WELLCOME RESEARCH LABORATORIES, TUCKAHOER, N. Y. 
Communicated August 29, 1944 


Sumner in the preceding article (pp. 285-294) has shown that the apparent 
concentration of guanine in the skin of Girella nigricans varies with the al- 
bedo of the background against which the fish are kept. 

Guanine was identified as a constituent of fish scales about eighty years , 
ago? and commercial “pearl essence’’ from this source is believed to con- 
sist essentially of free guanine.* The existence of guanine in fish scales may 
therefore be accepted, though with some reservations, since the chemical 
identification was carried out before isoguanine (2-hydroxy-6-aminopurine) 
which has properties nearly identical to those of guanine,* became avail- 
able. Within recent years a new series of pigments closely allied to the 
purines has been studied in detail by Wieland and associates.*® These 
include isoguanine (‘‘guanopterin’’)’ and a number of derivatives of the 
pteridine (azinepurine) heterocycle. These pigments were isolated first 
from butterfly wings, but at least one member now is known to exist in fish 
skin.’ Therefore, though the identification of guanine as a major constitu- 
ent of fish skin and scales may be accepted, the absence of similar sub- 
stances is by no meanis certain. 

The procedure adopted by Dr. Sumner for the determination of guanine, 
involves a preliminary separation of the material by a precipitation with 
copper and bisulfite, and finally a colorimetric estimation with phenol re- 
agent. The copper-bisulfite precipitation limits the number of substances 
which may be present to those of the purine group and substances closely 
allied chemically® but the reaction with the phenol reagent is relatively un- 
specific.’ Isoguanine, for example, and some of the pteridines might be 
expected to be determined as guanine by these procedures. Therefore a 
quantitative chemical or biochemical procedure was sought which would 
increase the specificity of the method. Ideally such a procedure, applied to 
small amounts of material, would distinguish between guanine and all 
other compounds known to behave like guanine with respect to copper 
precipitation and the phenol reagent. This would enable one to determine 
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what proportion of the material estimated as guanine was authentic gua- 
nine. Above all one could then determine whether the observed changes in 
concentration of the substance due to environmental conditions were due 
chiefly to changes in the concentration of guanine. To a considerable ex- 
tent a procedure including the use of the enzyme guanase’! was found to 
give the desired degree of specificity. This enzyme has been used by 
Schmidt’? for the determination of guanine in tissue extracts, but its action 
on substances closely related to guanine had not been studied previously. 
Consequently a study of the specificity of guanase action was undertaken. 

Guanase was prepared from acetone-dried rabbit liver powder by extrac- 
tion with water. The extract was purified by treatment with dialyzed iron 
and finally the enzyme was adsorbed on alumina cream,'* and eluated with 
a 0.2 M disodium phosphate solution. An amount of the substance to be 
tested equivalent to 0.1 mg. of amino nitrogen was incubated 18 hours at 
37° with guanase in 0.1 M phosphate buffer at pH 7.8. After the addition 
of an excess of sodium tetraborate solution, the ammonia which had been 


ACTION OF GUANASE ON THE PURINE OF Extracts oF Girella Nigricans 


GUANINE FOUND, MG. GUANINE PER G. FISH 


DIRECT ‘ AFTER COPPER PRECIPITATION 
SAMPLE* COLORIMETRIC COLORIMETRIC GUANASE 
Black 3.1 0.386 0.285 0.301 
3.2 0.347 0.250 0.245 
White 3.1 1.23 1.11 1.08 
3.2 1.04 0.91 0.85 


* Designations of Dr. Sumner (pp. 285-294). 


liberated was distilled into standard 0.02 N acid*and the excess acid was 
titrated. The action of guanase on twelve purines and azine-purines 
closely related to guanine was studied. Within the limits of experimental 
error, the enzyme showed no liberation of ammonia with the following 
substances: adenine, 7-methylguanine, 1, 7-dimethylguanine, 2-amino-6,8- 
dihydroxypurine (‘‘8-hydroxyguanine’’), 2-hydroxy-6-aminopurine (‘‘iso- 
guanine’), 2-amino-6-hydroxyazinepurine, xanthopterin, isoxanthopterin, 
xanthopterin carboxylic acid and isoxantopterin carboxylic acid. In 10 
experiments with guanine an average of 98 per cent of the theoretical 
amount of ammonia was liberated (range 93-106 per cent). Guanase 
showed a definite action on 1-methylguanine. In three experiments carried 
out as above an average of 81.0 per cent (maximum 84.2 per cent) of the 
theoretical nitrogen was liberated. The action of guanase was demonstrably 
slower on 1-methylguanine than on guanine itself. With a five-hour incuba- 
tion period and a somewhat feebler guanase solution, 69 per cent of the 
theoretical ammonia was liberated from guanine, while 8.5 per cent was 
being liberated from 1-methylguanine. 
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The action of guanase on the purine fraction of extracts of Girella ni- 
gricans then was studied. For this purpose fish were chosen to represent 
wide differences in apparent guanine content. The results of these experi- 
ments are given in table 1. Each extract contained about the same amount 
of material which reacted with the phenol reagent but was not precipitated 
by copper and bisulfite and thus was non-purine in character. In each in- 
stance the concentration of guanine found by the guanase method was 
identical with that found colorimetrically within the limits of accuracy of 
the enzymatic method. 

Whereas guanase failed to show a complete specificity in its action, by 
its use it was shown that isoguanine and a number of purines and pteri- 
dines could not be present in the extracts. The failure of guanase to dis- 
tinguish in all or none fashion between 1-methylguanine and guanine,‘ 
leaves open perhaps some possibility that the former substance might be 
present in the extracts, though 1-methylguanine has never been isolated 
from natural sources. To a considerable extent, this possibility was elimi- 
nated by a determination of the substance as the picrate.!® A solution of 
the fish purine containing 0.92 mg. total nitrogen, and 0.81 mg. guanine 
nitrogen as determined colorimetrically was precipitated at a volume of 6 
ml. by half saturation with sodium picrate. Titration of the precipitate 
gave an apparent guanine content of 0.85 mg. of nitrogen. When 1.00 mg. 
of guanine nitrogen was treated in the same way, 1.02 mg. were found by 
titration. On the other hand 1.00 mg. of nitrogen as 1-methylguanine gave 
a recovery of only 0.71 mg. due to the greater solubility of the picrate of 
this compound. 

The results leave little reason to doubt that substantially all the apparent 
guanine of the extract is authentic guanine, and that the amount of guanine 
in the fish skin changes with changes in the environmental conditions 
under which the fish are kept. 


* From The Wellcome Research Laboratories, Tuckahoe, New York. 
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germination of spores of Phycomyces blakesleanus, 1-methylguanine approaches in activ- 
ity guanine itself where other methyl guanines and closely related compounds have ac- 
tivities of a lower order of magnitude. (Unpublished experiments with Wm. J. Robbins.) 
Since methylation in the 1-position would result in fixation of the lactam tautomer, it is 
perhaps permissible to speculate that with guanase and with Phycomyces it is the lactam 
form of guanine which is active. 

8 Hitchings, G. H., and Fiske, C. H., Jour. Biol. Chem., 141, 827-835 (1941). 


ON SPONTANEOUS MUTATION 
By RICHARD GOLDSCHMIDT 
DEPARTMENT OF ZOOLOGY, UNIVERSITY OF CALIFORNIA 
Communicated August 3, 1944 


The occurrence of spontaneous mutation in a stock called px bl] was 
analyzed. Some of the facts have been reported previously,’ but there were 
some errors of interpretation in the preliminary report. The px Oi line is 
derived from standard plexus by mutation at the bs locus together with an 
appearance of modifiers for plexation and blistering. In an inbred and 
closely watched brother-sister culture an ‘“‘upheaval’’ (formerly called 
mass-mutation) occurred, in which mutants at the silver arc and rudi- 
mentary loci occurred simultaneously, together with return mutation of 
px and bs to wild. A similar upheaval was found twice later in mass cul- 
ture, the same loci being involved but in part with different alleles. All 
these mutants and return mutants appeared individually many times later 
in controlled broods of derivatives of px bl, or in in-and-out crosses of px bl 
derivatives. These mutants include a long series of characteristic alleles 
at the svr and a loci, repeated mutation at the bs and 0d loci, and further 
secondary return mutation from + to px and of the svr and a alleles to wild 
type. 

In a significant number of cases, two or more mutants or return mutants, 
including mutation from one allele to another, occurred simultaneously. 
This applies especially to the sur and a alleles. In such cases abnormal sex 
ratios were frequent. With the exception of these occurrences, the muta- 
tion rate was very low in px bi and its derivatives. But after in-or-out- 
crossing this rate rose considerably and also an accumulation of mutants 
was observed within a few generations after the first mutant had been 
found. This includes not only mutants and return mutants at the loci al- 
ready mentioned, but also the rather frequent occurrence of others at the 
Lobe, dachs, Notch, ebony loci. It is remarkable that in a somewhat small, 
tentative set of x-ray experiments an unusually high rate of mutation was 
observed, and among the mutants were alleles of sur and a as well as at the 
bs, N and d loci; of course, all lines were checked before the experiments 
for the absence of such mutants. 
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The px bi line as well as the lines of derived mutant stocks were thor- 
oughly analyzed genetically and cytologically. Aside from the presence 
of typical modifiers for the expression of the main loci, a series of non-recip- 
rocal, very small translocations from the X- into the autosomes was found, 
producing male lethal classes but otherwise viable. They did not contain 
known loci but could be localized approximately. The data show the inser- 
tion of small sections from near the left end of the X-chromosome into the 
second near arc and into the third near its right end. In addition, the ge- 
netic data require transpositions of the same X-chromosome loci, one into the 
right half between v-wy, the other to somewhere near the right end. Some 
but not all of these extremely small rearrangements could be found in the 
salivaries. The original px b/ stock already contained these features in dif- 
ferent combinations. There is a strong suspicion that a relation exists 
between the presence of these features and mutation at the sur and a loci. 
Strangely enough, some of the rearrangements were also found in an svr 
allele derived independently (in heat experiments) from Florida stock. A 
statistical check of the salivary chromosomes revealed the absence or the 
great rarity of conspicuous rearrangements for px b/ and the chief svr allele 
derived from px b/ in the first upheaval. But a very large number of 1-3 
band deficiencies and insertions could be found, some frequently, some 
seldom; they include the insertions from X- into autosomes already men- 
tioned. In the series of svr alleles the locus (left of the Bld break) attrib- 
uted to them through deficiency tests is always normal. In some of the 
alleles, however, a disturbance is typically found six bands to the right. 
In one allele it is a clear two-band inversion. In others the best possible in- 
terpretation of the abnormality is a one-band inversion. In the other 
long series of mutants at (or near) the a locus, only one shows a one-band 
deficiency. The other mutants derived from the line had either normal 
salivaries or were small deficiencies for known loci (Notch and black de- 
ficiencies). There is no reason to assume that the visible rearrangements 
were more frequent than the so-called point mutations. 

The facts in this closely followed case of spontaneous mutation, together 
with a detailed study of the lines involved, tend to show that mutation is 
not a haphazard phenomenon but the consequence of preéxisting conditions 
in the chromosomes which act upon other pointsin the same or other chromo- 
somes. Such facts—with which many data from existing literature may be 
considered as falling in line—suggest that mutation is not a chemical change 
within a gene-molecule, but a mechanical happening, probably a break 
under stress, produced by the synaptic properties of small rearrangements, 
including such breaks and rearrangements as are below the level of a single 
salivary band; i.e., so-called point mutants. The latter can now be en- 
visaged after the intimate structure of one band has been elucidated by 
Kodani.’ A full report is ready for press. 
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1 Spontaneous chromatin arrangements and the theory of the gene, Proc. Nat. Acad. 
Sct., 23, 621-623 (1932). 

2 Kodani, M., ‘“‘The Structure of Salivary Gland Chromosomes of Drosophila melano- 
gaster,”’ Jour. Hered., 33, 114-133 (1942). Goldschmidt, R., and Kodani, M., ‘‘The 
Structure of the Salivary Gland Chromosomes and Its Meaning,” Am. Nat., 76, 529-551 
(1942). 


ONE-BAND INVERSION 
By RICHARD GOLDSCHMIDT AND ALOHA HANNAH 
DEPARTMENT OF ZOOLOGY, UNIVERSITY OF CALIFORNIA 


Communicated August 3, 1944 


In a recent paper by Oliver and Green! the following remarks are found: 
“In some of the females heterozygous for /z*, some but not all of the cells 
of the salivary glands show a torsion of one heavy band which is suggestive 
of a small inversion or possibly a duplication.’”’ The feature was found only 
in the presence of the allele /z*. It happens that we have been studying in 
detail what probably is a parallel case. The mutants in question are a 
series of alleles at the silver locus (1-0.1) of Drosophila melanogaster which 
were studied closely in connection with work on spontaneous mutation. 

On the basis of deficiency tests the silver locus is expected to be located 
to the left of 1C in the first chromosome. In all alleles, including those pro- 
duced by x-rays, this region is perfectly normal. But in some alleles the 
region farther to the right, between the B/d break and the bulb, shows typi- 
cal disturbances. The only one of these which can be described with cer- 
tainty is found in the allele svr?’*. Here a 2- (4-) band inversion of 
1E1-4 is always found. These bands are described in Bridges’ standard 
map (1938) as two double bands, the one to the left being thinner, that to 
the right being rather thick. We think that there exist actually only two 
bands, the double nature of which is not always visible, being based upon a 
partial separation of the two groups of perultimate chromomeres of which 
the bands are constituted (see Kodani, 1942; Goldschmidt and Kodani, 
1943).2. In the heterozygote this inversion looks as pictured in figure 1. 
One cannot always be sure that the interpretation as a minute inversion is 
correct. But in the homozygote the order of the thin and thick bands 1-2, 
3-4 is clearly reversed, as shown in figure 2. Therefore we feel confident 
that the interpretation is correct. 

In some of the other alleles, namely, sur”, sur?% *%, sur? 4h, an ab- 
normality of the same bands is found in a majority of X-chromosomes (this 
is never true of the controls), a structure which offers great difficulties of 
interpretation. The order of the bands 1-2, 3-4 is normal. But the thick 
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band 3-4 is distorted in a way which is difficult to describe. When focusing 
high and low the band seems normal, but in between a kind of cross or figure- 
eight appears, an apparent crossing of 1, 2 with 3,4. This seems to be the 
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result of two features; first, the arrangement of the chromatin on the sur- 
face of the disc-like band (like a ring), and secondly, a dislocation of the 
two halves of this ring belonging to the two homologous chromosomes, 
by a twist in opposite directions. This means that viewed from 
above the upper rim of the disc in one chromosome is shifted to the right, 
and that of the other chromosome to the left, so that a kind of figure-8 
becomes visible. This is pictured in figures 3 and 4, where the same chro- 
mosome is drawn at different foci. Our interpretation of these typical as- 
pects, which can hardly be drawn as they appear in focusing, is that the 
distortion within the disc represents a one-band inversion. This requires 
the assumption that the left and right surfaces of a single disc are serially 
different. The new discoveries on the structure of the individual disc (see 
Kodani, and Goldschmidt and Kodani*) make such an interpretation pos- 
sible, since they revealed the presence of a coil of genonema with 
an attached series of perultimate chromomeres in a disc (each one double). 
An inversion of one disc would change the direction of the coil, creating a 
tension in the state of synapsis. (Consult figure 6 in Kodani’s paper.) 
Homozygous chromosomes were always normal. 

In sor * a different structure is found, again for the same bands. We 
did not succeed in finding a satisfactory interpretation. The only certain 
fact is that these bands are not normal. Drawings will be given in a forth- 
coming full report. 

If our interpretation is correct—and we have no doubt that it is—the 
group of mutants which turn out to be small rearrangements will probably 
be considerably enlarged. Simultaneously the senior author’s interpreta- 
tion of the experimental work on the structure of the salivary chromosomes? 
will receive important support. 


1 Oliver, C. P., and Green, M. M., ‘‘Heterosis in Compounds of Alleles in Drosophila 
melanogaster,’’ Genet., 28, 331-387 (1944). 

* Kodani, M., ‘““The Structure of Salivary Gland Chromosomes of Drosophila melano- 
gaster,’’ Hered., 33, 115-133 (1942). Goldschmidt, R., and Kodani, M. The Structure 
of the Salivary Gland Chromosomes and Its Meaning, Amer. Natur., 76, 529-551 (1942). 
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CHROMOSOME ASSOCIATIONS IN TETRAPLOID HYBRIDS 
BETWEEN PRIME TYPE 1 AND PRIME TYPE 2 IN DATURA 
STRA MONIUM 


By A. D. BERGNER* 


UNITED STATES DEPARTMENT OF AGRICULTURE, SPECIAL GUAYULE RESEARCH PROJECT, 
SALINAS, CALIF. 


Communicated September 8, 1944 


Prime types among diploid plants of Datura stramonium have been dis- 
cussed in several previous publications.':*:* They are the result of inter- 
changes between non-homologous chromosomes. No studies have been 
made of chromosomal configurations in the F, hybrids between tetraploid 
prime types (PT’s) in Datura with the exception of PT 1 X PT 2. 

In diploid races which are homozygous for PT’s, the 24 chromosomes are 
associated as 12 bivalents at the meiotic first metaphase (/ I) in pollen- 
mother-cells (PMC). If homozygous PT 1 (standard Line 1) is crossed 
with any other homozygous PT, the chromosomes of the hybrid show a 
configuration of four or more chromosomes plus 10 or fewer bivalents, de- 
pending upon which PT is crossed with PT 1. In tetraploid races which 
are homozygous for PT”s, the 48 chromosomes are associated as 12 quadri- 
valents, or the equivalent, at M I since there are four of each kind instead 
of the two in diploids. In a tetraploid hybrid between PT 1 and any other 
PT many different types of configuration are possible. 

Diploid Prime Type 2.—PT 2 has the chromosomes 1-184 and 2-17 in- 
stead of the 1-2 and 17-18 chromosomes of PT 1. The hybrid between 
PT 1 and PT 2 shows a © 4 + 10 bivalents (BV) at MI. The chromo- 
somes in the © are as follows: (1-2-2-17-17-18-18-1).6 Occasionally 
instead of the © 4, which is illustrated in figure A, these four chromosomes 
form a chain (ch) of -18-1-1-2-2-17-17-18 or more rarely -2-1-1-18-18-17- 
17-2. They may also occur as the two pairs -2-17-17-18 + 18-1-1-2. 
The © 4 and ch 4 are usually so oriented on the spindle that disjunction is 
alternate. 

Under favorable conditions these arrangements can be determined be- 
cause of size differences among these four chromosomes. The L-sized 1-2 
chromosome is twice as long as the m-sized 17-18 chromosome; 1-18 
is slightly longer than 17-18 and 2-17 is slightly shorter than 1-2. These 
differences in length were first determined by using the location of the extra 
chromosome in certain 2n + 1 types heterozygous for PT 2 at M I in 
PMC’ and later checked by measuring these chromosomes in pollen grains.® 
The 1-18 and 17-18 chromosomes probably have the same insertion region; 
correspondingly the 1-2 and 2-17 chromosomes have the same insertion 
region. This conclusion is based on a comparison of total length and length 
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of each arm of these four chromosomes in PMC and in pollen grains. The 
-18 arm of the 1-18 chromosome, therefore, is the same as the -18 arm of 
the 17-18 chromosome and the -2 arm of the 2-17 chromosome is the same 
as the -2 arm of the 1-2 chromosome. Concerning the -1 arm of the 1-18 
and the -17 arm of the 2-17 chromosomes, assumptions have been made 
which are based on the appearance of the 2m + 1-18 and 2m + 2-17 ter- 
tiary plants. The -1 arm of the 1-18 chromosome consists of approxi- 
* mately the distal two-thirds of the -1 arm of the 1-2 chromosome. The 
-17 arm of the 2-17 chromosome consists, distally, of the -17 armofthe17-18 
chromosome plus, proximally, of the inner third of the - 1 arm adjacent to the 
insertion region. This is diagrammed in figure 1 where the 1-2 chromosome 
is represented by the solid line, the 17-18 chromosome by the dashed lined. 

Tetraploid Prime Type 2.—In 1930’ it was 
reported that a © 8 wasseeninthe4n (PT 1 . 2 
1 X PT 2) hybrid. At the Seventh Inter- 
national Genetics Congress’ the © 8 and 
some other types of configuration were 
shown. It was suggested by Dr. F. W. 
Sansome that the relative frequency of the , n 18 
different possible types of configuration. 
should be studied. Before this could be FIGURE 1 
done, however, a new supply of 4n homo- Lines representing approxi- 
zygous PT 2 plants was needed. By soaking mate portions of the 1-2 and 
seeds from 2n homozygous PT 2 plants in a 1718 chromosomes of PT 1 in 
: rete the 1:18 and 2:17 chromosomes 
solution of colchicine, such 4n plants were of PT 2. Solid line denotes 1-2 
obtained. From the cross of 4n PT 1 X 4m chromosome, dashed line the 
PT 2 there were obtainéd hybrids which 17-18 chromosome. 
were grown to maturity. The present report 
is based on observations made on 45 PMC from ten of these plants. 

Among the 48 chromosomes of 4n PT 1 there are four 1-2 and four 17-18 
chromosomes; likewise for 4n PT 2 there are four 1-18 and four 2-17 
chromosomes. In the 4n (PT 1 X PT 2) hybrid there are two each of the 
1-2, 2-17, 1-18 and 17-18 chromosomes. Whenever it is possible to iden- 
tify these chromosomes, their arrangement at M I can be determined. 
Drawings of three of the many possible types of configuration involving 
these eight chromosomes are shown in figure 2, B—D, In B the sizes of the 
chromosomes in the © 8 are indicated. In C the four L-sized chromosomes 
form a © 4 to which is attached a ch 4 consisting of m-sized chromosomes. 
This configuration can be diagrammed as follows: (1-2—2-17)s"°-17-18- 
18-1-1-18-18-17. In D there are two ©’s 4, each of which is (1-2-2-17- 
17-18-18-1). 

In table 1 are listed the chromosomal arrangements seen in the 45 PMC 
referred to above. This number is admittedly too small for a study of rela- 
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TABLE 1 
Types oF CONFIGURATION INVOLVING THE 1-2, 1-18, 2-17 AND 17-18 CHROMOSOMES AT MEIOTIC 
First METAPHASE IN 45 POLLEN-MOTHER-CELLS FROM TEN TETRAPLOID F; Hysrips 
(PRIME TYPE 1 X Prime Type 2) oF Datura stramonium 


CHROMOSOMES IDENTIFIED AS 
(1-2)2 + (2-17)2 + (17-18)2 + (1-18)2 
(1-2 — 2-17). + (1-18). + (17-18)2 


(2-17 — 17-18). + (1-2)2 + (1-18)2 
(1-18 — 18-17). + (1-2)2 + (2°17)2 


(1-2—2-17—17-18—18-1)+(1-2—2-17—17-18—18-1) 
(1-2—2-17—17-18—18-1)+(1-2—2-17—17-18—18-1) 
(1-2—2-17—17-18—18-1)+(1-2—2-17—17-18—18-1) 


(1-2 — 2-17). + (1-18 — 18-17). 


(1-2 — 2-17). + 1-18 — 18-17 — 17-18 — 18-1 


(2-17 — 17-18 — 18-1)2 + (1-2)2 
(18-1 — 1-2 — 2-17). + (17-18)2 
Chain of 6 + (17-18)2 

1-2— 


1-2 — (2-17). — 17-18 — 18-1 + 1-18 — 18-17 


Sequence not recorded 
Sequence not recorded 
Sequence not recorded 
Sequence not recorded 


Sequence not recorded 


2-1—1-18—18-17—17-18—18-1—1-2—2-17—17-2 


(1-2). — (2-17 — 17-18 — 18-1): 
(1-2), — (2-17 — 17-18 — 18-1), 
(2-1), — (1-18 — 18-17 — 17-2) 


(1-2). — 2-17 — 17-18 — 18-1 — 1-18 — 18-17 — 17-2 
(2-1). — 1-18 — 18-17 — 17-2 — 2-17 — 17-18 — 18-1 
(2-17): — 17-18 — 18-1 — 1-2 — 2-1 — 1-18 — 18-17 
(17-18), — 18-1 — 1-2 — 2-17 — 17-2 —2-1— 1-18 


(1+2)2 — (2-17 — 17-18 — 18-1)s 


2-17 — 18-1 
(1-2), — 2-17 — (17-18). — 18-1 
Unresolved configurations of 8 
Unresolved configurations of 8 
Unresolved configurations of 8 


@ THEORBTICAL NUMBER OF ARRANGEMENTS. 
6b ASSOCIATIONS OF THE OTHER CHROMOSOMES. 
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tive frequencies of the many possible types of configuration. In fact, not 
all of the possible types were found. Nor do these 45 PMC represent a 
random sample because no PMC was included if one or more univalents 
were present. This extreme precaution was taken to make certain that 
the plants were true tetraploids, having four each of the following ten kinds 
of chromosomes: 3-4, 5-6, 7-8, 9-10, 11-12, 13-14, 15-16, 19-20, 21-22, 
23-24 and two each of the 1-2, 1-18, 2-17 and 17-18 chromosomes. 

In column 1 of table 1 are listed in diagrammatic form thirteen types of 
configuration found in 39 PMC. In column 2 the theoretically possible 
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FIGURE 2 
Chromosome associations at meiotic first meta- 
phase in pollen-mother-cells of 2m and 4n (PT 
1 X PT 2) plants of Datura stramonium. A. 
2n: © 4 (LLmm) + 10 BV. B. 4n: ©8 
(LmmLLmmL) + 10 quadrivalents. C. 4n: © 
4 (LLLL) with attached ch 4 —mmmm + 9 
quadrivalents + 2 BV. D. 4n: two ©’s 4 
(LLmm); other chromosomes omitted. 


number of arrangements in each of these types is given. Whenever the 
1-2, 1-18, 2-17 and 17-18 chromosomes could be identified, as was pos- 
sible in 26 PMC, the sequences of the chromosomes in the configurations 
could be determined. They are listed in column 3. Further separations 
were made according to the number of quadrivalents and bivalents formed 
by the 40 other chromosomes. They are indicated in column 5 where the 
subscript ;y denotes quadrivalents and the subscript ;; denotes BV. These 
latter data are included to show to what extent quadrivalents were formed 
in the same PMC under consideration. This tendency to form quadriva- 
lents is consistent with the known behavior of other 4m Daturas. 
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Of the thirteen types listed, nine are distinct. Numbers 4, 6, 9 and 11 are 
modifications of Nos. 3, 5, 8 and 10, respectively. They form a series ex- 
tending from four BV through a © 4 + 2 BV, two ©’s4, ©6+ BV toa 
© 8. One quadruple attachment transforms the © 8 toa © 6 attached to 
a BV. 

Discussion.—In the 2n (PT 1 X PT 2) hybrid, since there is only one 
each of 1-2, 1-18, 2-17 and 17-18, a © 4 or ch 4 is the largest possible con- 
figuration and it is found regularly except when there are two pairs. In 
the 4n (PT 1 X PT 2) hybrid, however, there are two each of these chro- 
mosomes. Consequently one may find four BV. In addition many types 
of configurations which involve four, six or eight chromosomes are possible. 
Even in so modest a sample as 39 PMC, nine types plus four modifications 
were found. 

Table 1 shows that practically two-thirds of the configurations involved 
all eight chromosomes. Moreover the © 8 was the type found most often 
(24 per cent). By including the ch 8 also the proportion was 31 per cent. 
Seven PMC (16 per cent) showed a © 6 or ch 6 attached toa BV. These 
were the commonest alternatives to a © 8 and involved one multiple at- 
tachment. Closely related is the © 6 or ch 6 plus a BV, where the BV is 
not attached to the other six chromosomes. If Nos. 5, 6, 10 and 11 are 
combined, their total equals 22 per cent. 

Configurations involving four chromosomes were second in frequency. 
Two ©’s4each or a © 4 + ch 4 were found in 8 PMC or almost 18 per cent. 
One © 4 + 2 BV amounted to 7 per cent. Combining these two types 
gives a total of over 24 per cent. 

Most of the configurations which involved four ur more chromosomes 
were merely concatenations since they consisted solely of a series of at- 
tachments at the homologous ends of only two chromosomes. Theoretically 
three or four chromosomes could be attached together at their homologous 
ends. In 10 PMC there were 14 instances of such multiple attachments. 
Six triple attachments are listed under Nos. 7, 11 and 13 in table 1. Eight 
quadruple attachments are listed under Nos. 7, 10, 12 and 13. Three PMC 
showed more than one multiple attachment: No. 7 showed one triple and 
one quadruple, No. 12 showed two quadruple and No. 13 showed one triple 
and two quadruple attachments. 

Seven PMC are classified under Nos. 10 and 11 in table 1. The multiple 
attachment involved a BV to which was attached a © 6 or ch 6. In five 
_ PMC the 1-2 chromosomes formed the BV, in one the 2-17. Also in No. 12 
the 1-2 chromosomes formed a BV. Possibly because of their large size 
the 1-2 and 2-17 chromosomes may form a BV more frequently than the 
m-sized ones. 

It is to be regretted that the sequences of the chromosomes in the many 
©’s 8 were not recorded with sufficient certainty to warrant publication. 
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In one of the © 8 + 10;y PMC the sequence probably was (1-2-2-1- 
1-18-18-17-17-2-2-17-17-18-18-1). In one of the © 8 + 9:y + 29, PMC 
the sequence probably was (2-1-1-18-18-17-17-18-18-1-1-2-2-17-17-2). 
This same sequence is given for one of the ch 8 + 10;y PMC. 

In the 2” hybrid only one kind of © 4 can be formed, namely (1-2-2-17- 
17-18-18-1). In the 4m hybrid, however, in addition to this arrangement, 
other arrangements in a © 4 are possible because each of these chromo- 
somes is present in duplicate. It is possible to have a © 4 which consists 
of two of each of two kinds of chromosomes. These possibilities are as 
follows: (1-2-2-17)2, (1-18-18-17)2, (2-17-17-18)2 and (2-1-1-18):. This 
last kind was not found at all and (2-17-17-18)2 was found only once in a 
© 4+ 2 BV type of configuration. 

There are three possible arrangements for two configurations of four 
each as follows: (2-1-1-18-18-17-17-2) + (2-1-1-18-18-17-17-2), (1-2- 
2-17). + (1-18-18-17), and (2-1-1-18), + (2-17-17-18)2. In seven PMC 
two ©’s 4 were found and in another a © 4 + ch 4. If these two types are 
considered together, it is seen that the first two arrangements were identi- 
fied and found with equal frequency. However, the third possible arrange- 
ment was not identified. 

Summary.—In tetraploid Datura hybrids.between prime type 1 and prime 
type 2 there are two each of the 1-2 and 17-18 chromosomes from PT 1 
and of the 1-18 and 2-17 chromosomes from PT 2 and four each of the 
other chromosomes. 

In 45 PMC from ten plants the configurations formed by associations 
among the 1-2, 1-18, 2-17 and 17-18 chromosomes at M I were studied. 
Two-thirds of the configurations involved all eight chromosomes. 

In 39 PMC the configurations could be diagrammed. Nine distinct 
types plus four modifications were found. The series extended from four 
BV through © 4+ 2 BV, two ©’s4,06+ BVtoa©8. In 24 per cent 
of the PMC the © 8 was found. A common modification was a ch 8. 
A © 6 or ch 6 attached to a BV also was seen frequently. Two configura- 
tions of four each were the next commonest types. 

Usually a configuration which involved four or more chromosomes con- 
sisted solely of a series of attachments at the homologous ends of only two 
chromosoms (concatenation). Only fourteen instances of attachment at 
the homolo sous ends of three or four chromosomes were found. 

The 1-2 chromosomes formed the BV in five of the seven PMC in which 
a © 6or ch 6 was attached toa BV. 

In 2n hybrids between PT 1 and PT 2 only one kind of chromosome ar- 
rangement is possible in the © 4 (1-2-2-17-17-18-18-1). In 4n hybrids 
three additional arrangements in a © 4 were found, namely (1-2-2-17)s, 
(1-18-18-17)2 and (2-17-17-18)s. 
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* This work was done at the Carnegie Institution of Washington, Department of 
Genetics, Cold Spring Harbor, New York, while the author was associated with Dr. A. 
F. Blakeslee. 
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ON MECHANICAL SELF-EX CITED OSCILLATIONS 
By N. MINorsky 
Davip TAYLOR Mopet Basin, Navy DEPARTMENT 


Communicated August 28, 1944 


1. Introduction.—Self-excited oscillations, particularly electrical ones, 
form the object of numerous studies in recent years. Mechanical oscilla- 
tions of this kind have been less explored although they are frequently ob- 
served in practice. 

The object of this note is to describe a typical phenomenon of this kind 
observed during experimental work on the antirolling stabilization of 
ships by the activated tanks method. The problem of the antirolling sta- 
bilization of ships itself does not form the object of the note and only a few 
words will be sufficient to explain the principal details of these experiments. 

2. Experimental Arrangement and Principal Data.—The principal fea- 
tures of ship stabilization by tanks can be conveniently studied both theo- 
retically and experimentally by means of two coupled pendula of which 
one P is a rigid physical pendulum and the other P’ is a pendulum consist- 
ing of a liquid (water, oil, etc.) filling a U-tube connected to two tanks and 
secured to the rigid pendulum P so as to obtain an oscillation coplanar with 
the latter. : 

In the experiments described here the activation of the ballast was ac- 
complished by means of an axial variable pitch pump inserted in the U- 
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channel of the liquid pendulum. The impeller pump was kept running at 
a constant speed and the activation was accomplished by controlling the 
blade angle a of*the runner by a control equipment responsive to the 
angular motion of P. 

Let @ be the angle of P and ¢ the relative angle formed by the levels of 
P’ with respect to P. The rate of flow of the ballast is then proportional to 
¢. The following experimental results will be of importance. 

1. It was found that the roll-quenching efficiency is greatest when the 
following condition of control is fulfilled 


a = 6 (1) 


that is, when the blade angle a is made to vary in proportion to the angular 
acceleration 6 of the rigid pendulum. The coefficient of proportionality ) 
could be adjusted by changing the amplification factor of an electronic cir- 
cuit controlling the blade angle a. 

2. For a steady harmonic actuation of the ballast by the impeller pump 
it is observed that there exists the following relation: 


g= 6 (2) 


between the rate of flow ¢ and the angular acceleration 6; this means in 
view of (1) that the rate of flow ¢ and the blade angle a are approximately 
proportional to each other or, 


a=h/e = ¥ (3) 


the coefficient e changes somewhat according to conditions. 

3. Moreover, from numerous tests, it has been ascertained that the 
force imparted by the impeller pump to the liquid column increases more or 
less linearly with a for small values of a and increases less than in propor- 
tion for larger values of a. The moment M of the force created by the 
impeller and applied to P’ can be, therefore, approximated by a relation 


M = Mia — Mga’ (4) 


where M/, and Ms; are independent of a. The form of the curves approxi- 
mated by this expression varies somewhat for different conditions. 

3. Phenomenon of Self Excitation —It has been ascertained that when 
the coefficient \ of amplification is made large enough, a self-excited os- 
cillation of relatively high frequency sets in and superimposes itself on the 
normal stabilizing action. If the coefficient \ continues to increase further, 
self-excited oscillation exhibits the presence of harmonics. If, however, 
d is decreased gradually, at a certain point \ = Xp the self-excited oscilla- 
tion disappears. During the experiments the parasitic fluttering of the 
blades had a frequency nearly five times that of theoscillatory process 
which the equipment is intended to quench. 
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The interesting feature of the phenomenon lies in the fact that it may 
originate also when the pendulum P is almost at rest. In fact it is suf- 
ficient to disturb the pendulum at rest very little, e.g., by touching it 
slightly, to be able to release violent surges of water in the tanks with the 
corresponding blade angle oscillations, which persist indefinitely. The 
pendulum P also begins to oscillate with the same high frequency but 
with an amplitude of a fraction of one degree whereas the oscillation of P’ 
is generally very appreciable and with strong amplification may reach the 
value of about 20 to 25 degrees. In what follows we shall investigate the 
problem in this particular case which is relatively simple. 

4. Differential Equations——The fact that during the self excitation 
from rest the amplitude of motion in the 6 coérdinate is very small, makes it 
possible to neglect the effect of the passive component of coupling between 
P and P’ since this component depends on the existence of a finite oscilla- 
tion in the 6 degree of freedom, and consider only the active component of 
coupling, i.e., the action of the moment VM due to the pump on P’ and the 
reaction of this moment on P. 

The simplified differential equations of the system in this case are: 


where j, k and c for the liquid pendulum P’ and J, K and C for the rigid 
pendulum P are certain constants, the physical significance of which is suf- 
ficiently obvious. The left-hand terms of these equations are linearized 
somewhat which has no further influence on what follows. The neglect of 
the passive component of the coupling as just mentioned is in the fact that 
the first equation contains only the relative codrdinate ¢ instead of the ab- 
solute one g — @ since in this particular case @ is very small. 

Equations (5) are thus written down with the knowledge of the principal 
features of the phenomenon and our purpose will be to fit these features into 
the form of solutions we are going to obtain. Using the equations (4) and 
(3) we obtain 


je + (k — m)¢ + my* + co = 0 (6) 
I6+ K6+ Co = —mg + me (7) 
where m, = M;d/e and ms; = M3\3/eé*. (8) 


It is seen that owing to the neglect of the passive component of coupling 
equations (6) and (7) do not form actually a coupled system since equation 
(6) contains only y. Hence, if we succeed in solving the non-linear equa- 
tion, ¢ will appear as a known function in (7) which will permit determining 
6 and through the equation (1) of control we shall be able to establish the 
condition of coupling. 

Thus, the principal interest centers on the equation (6). It is noted 
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that for a sufficiently small value of \, the coefficient m, is also small so 
that k — m,> 0. In such a case it is apparent that no self excitation is 
possible in the ¢ degree of freedom since equation (6) shows that the liquid 
pendulum P’ has a purely. dissipative damping of a combined velocity, and 
velocity-cube type. The really interesting case arises when the parameter 
\ is increased so as to make the coefficient k — m, < 0 which case we shall 
consider from now on. Physically this means that the energy imparted 
by the blades outweighs the dissipation of energy caused by the velocity 
damping k¢. 

In order to bring the problem within the scope of non-linear mechanics 
we shall assume that the damping terms are relatively small in comparison 
with the inertial (j¢) and stability (cy) terms which is sufficiently justified 
by the actual order of magnitude of these quantities. This permits writing 


m —k = ups; Ms; = uq 


where y is a small parameter (u > 0) and s > 0 and g> 0 are certain con- 
stants of the order of j and c. With these conventions equation (6) can be 
written : 


je + ce + u(—se + ge") = 0. (9) 


This equation is of the Rayleigh type. By differentiating it becomes of the 
Van der Pol type in ¢ and thus may be shown! to have the approximate 
solution 


1 che 
g=— splitece COS (wot + Yo) 


Wo Ms 


g=- (ea sin (wot + yo). (10) 

Substituting this value for ¢ into equation (7) one easily obtains @ since 
this equation is linear. Since the term m;¢* contains the third harmonic it 
is apparent that the oscillation of P will also contain that harmonic and in 
view of equations (1) and (3) this harmonic will also appear in the ¢ oscilla- 
tion. The amplitudes of these harmonics are small if the ratio m3/m, is 
small as we have assumed. We shall limit this study to the first approxima- 
tion only assuming that m; = 0 in equation (7). Since the phase y is 
clearly of no interest here we can write 


Il6+ K6+ CO = —meo =m a) sets —® ai = Ae 
3m3 . 


4(m, = k) 


where A = edt 3 ; whence the approximate expression for the 
Ms 
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amplitude 6 is 
at Tex? 3mz3 


and the amplitude of the angular acceleration is then 


; _ m on — k) 
see ae: es ay) 


The blade angle is given by equation (1) 
a = rb. 


It is seen that the blade angle generally increases with the amplification 
factor but not exactly in proportion to A since by equations (8) both m 
and ms; are functions of \. It is generally observed that the blade angle a 
increases with the amplification factor but the observation is handicapped 
owing to the presence of harmonics which are left out in this study. 

6. Critical Value of Parameter—We shall now attempt to establish 
analytically the principal feature of the phenomenon mentioned at the be- 
ginning of Section 3, namely, a rather sudden appearance of self excitation 
at a certain value X = Xo of the parameter when it is increased from small 
values as well as its sudden disappearance at the same value \ = Xo, when A 
is gradually decreased from larger values. 

Although this fact has been discussed to some extent in Section 4 we shall 
analyze it now from a more general standpoint resulting from the theory 
of Poincaré-Liapounoff.? Only equation (6) will interest us here because 
equation (7) does not present any difficulty. 

Putting m, = da, and ms; = \*a3 where a; and a; are certain positive con- 
stants at least within the range of the experiment as follows from (8) we 
write equation (6) in the form 


je + (& — dane + mag? + cp = 0. (6a) 


This equation is equivalent to the following system of differential equations 


ey 
y = —g(\)y — woe — ms’y? (12) 


where 


g(a) = (k — Aay)/j and m;’ = m;/j. 


We shall investigate the conditions of stability of the liquid pendulum 
when at rest. Liapounoff has shown that for this purpose it is legitimate 
to neglect the non-linear term in the second equation (12) which gives 
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ede 
y = —g(r)y — we. (13) 


We can now apply the classical method of Poincaré which consists in in- 
vestigating the singularities of the system (13) as a function of the pa- 
rameter. In this case there is clearly one singular point of interest, viz., 
¢g =y=0. The characteristic equation is 


S? + g(A)S + wo? = 0 (14) 


and its roots are 
2 
Si,2 = —g0)/2 * em — a" (15) 


There are four possibilities according to the sign of g(A)/2 and of the 

quantity under the square root. 
| | 

If ee > wo the roots are real and of the same sign which corresponds 
to a nodal point of the system. This case is to be discarded because the 
experiment always shows that a self-excited oscillation builds itself up in a 
series of swings and never reaches a steady state aperiodically which would 
require a pump of an enormous power. Hence, we have to consider the 


case when os < wo in which case the roots are conjugate complex and 





the singularity ¢ = y = 0 is a focal point, stable if g(A) > 0 and unstable if 
g(A) < 0. The value 4 = Ao for which g(A9) = 0 and which separates the 
stable focal points from the unstable ones is thus the critical or bifurcation 
value of the parameter \ (Poincaré). 

Since g(A) is monotone decreasing with A increasing and is positive for 
small ), it follows that by gradually increasing the value of the parameter 
\ the original stable focal point for \ < A» becomes unstable for \ 2 do. 

From the theory of Poincaré it follows that for the critical value \ = Xo 
for which a stable singularity becomes unstable there may appear at least 
one stable limit cycle which means a stationary oscillation. The proof of the 
existence of a limit cycle is rather laborious and requires the investigation of 
a series expansion satisfying the non-linear system (12). This need not be 
done here since we know that the self-excited oscillation exists for \ > Yo 
and the theory of the first approximation gives the answer in a simple 
manner. 

7. Conclusion.—It is seen that the principal features of the phenomenon 
in question are consistently explained on the basis of non-linear mechanics. 
In fact, if one had to follow purely linear methods it would be impossible 
even to predict the existence of a phenomenon of this nature. 
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FIRST PROOF THAT THE MERSENNE NUMBER 
My; IS COMPOSITE 


By H. S. UHLER 
SLOANE Puysics LABORATORY, YALE UNIVERSITY 
Communicated September 9, 1944 


The form of a Mersenne number is M, = 2? ~- 1 where p is prime. Just 
three hundred years ago Mersenne published, in effect, the statement that 
the only values of p not greater than 257 which make M, prime are 2, 3, 5, 
7, 13, 17, 19, 31, 67, 127 and 257. The prime or composite character of all 
of the 55 numbers included under this conjecture, except the six corre- 
sponding, respectively, to = 157, 167, 193, 199, 227 and 229, had been in- 
vestigated prior to the year 1935. Contrary to Mersenne’s surmise it was 
found that for p = 67 and 257 M, is composite, and that M, is prime for 
b = 61, 89 and 107. The data presented above have been derived from a 
comprehensive paper by R. C. Archibald.! 

The present contribution marks the first fruits of a friendly suggestion 
made last year by Professor Archibald and seconded by Professor D. H. 
Lehmer that the author turn his attention to problems of factorization in 
general, beginning with the special problem of the character of the six 
heretofore uninvestigated Mersenne numbers. 

The modern technique of this problem is based explicitly upon the fol- 
lowing theorem, discovered by E. Lucas and clarified by D. H. Lehmer,? 
namely: “The number N = 2" — 1, where n is an odd prime, is a prime if, 
and only if, N divides the (n — 1)-st term of the series, 


= 4, Ss = 14, AR ab ,- Sage rae 


where S, = Sp—1 — 2.” 

Since the number Mj57, which is the subject of the present study, equals 
182 68770 46663 62864 77546 06040 89535 37745 69915 67871 and since the 
above theorem requires that this 48-digit number be used as a divisor 149 
times (8 to 156, inclusive) it should be obvious that the prospective inves- 
tigator would focus attention upon this aspect of the computation. The 
procedure followed by R. E. Powers® in showing that Mo is composite ap- 
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pealed at first to the author and hence it was used as far as the 34th member 
of the Lucasian sequence. At this stage Powers’ technique was abandoned 
because it had proved to involve an inordinate amount of writing and men- 
tal arithmetic and especially because it did not make use of a computing 
machine exclusively. 

The patent idea of substituting for division multiplication by the re- 
ciprocal of the divisor overcame all of the objections entertained by the 
writer against all of the schemes for direct machine-division which he had 
read in the literature. Consider s,’ = N-q + 7r;. Then (1/N)-s;’ = & 
+ r;,/N. Hence if we multiply any term (>) of a congruently equivalent 
Lucasian sequence by 1/N the succession of figures to the left of the decimal 
point will represent the integer g,, while the digits to the right of the deci- 
mal point will constitute an approximation to the incommensurable decimal 
fraction 7,;/N. Theoretically 7, may be obtained by multiplying r,/N by 
N without evaluating a single g, at any stage of the work. This procedure 
was not followed because it involves a greater number of figures than the 
alternative process of multiplying g, by N and then subtracting N-q, from 
s,’, that is, defining r, as 5,’ — N-q,. 

The following details of the actual computation should make the matter 
perfectly clear. Ribbon A of cross-ruled paper has transcribed in nonads 
along the lower edge a sufficient number of figures of the reciprocal of 
N(Mis7) thus: 1/N = 0-(47 zeros) 547382212 626881668 329581868 
472623488 781352321 048510246 ---. Similarly the bottom line of ribbon B 
reads: N = 182687704 666362864 775460604 089535377 456991567 
871.000000. Ribbons A and B do not have to be changed throughout the 
work from k = 8 tok = 156. On the other hand, each value of k requires 
three ribbons of its own. Along the top line of a longer ribbon of cross- 
ruled paper is written, for example, the value of 755 as 744260418 487541370 
275270351 147668949 680082908 44.0000000. This number was machine- 
squared by D. N. Lehmer’s‘ algorithm and recorded along the lowest line 
of the same ribbon. After subtracting 2 this line read sis.’ = 553923570 
527250212 304427560 946230012 742643893 708358589 871125964 
130378182 986908821 472209423 2334. The undiminished squared num- 
ber was carefully checked with the moduli 10* + 1, 104+ 1 and 10’ + 1. 
Next ribbon A was laid on the longer ribbon so that the nonads of 1/N 
were juxtaposed above and in exact register with the first six nonads of 5156’. 
The multiplication of sys6’ by 1/N was then effected and the resulting value 
of the integer giss presented itself on the bottom line of a fourth strip of 
paper. This short act of multiplying was always repeated as a preliminary 
check on its accuracy. Actually gu = 303207909 661388759 507868593 
612664030 922503212 07. Finally ribbon B was aligned with the nonads 
of diss and a portion of the product N-gis¢ was recorded on the second line 
from the bottom of a fifth strip of cross-ruled paper. This multiplication 
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was performed only on the nonads from the fifth to the eleventh, inclusive, 
since the first four sums-of-products would have reproduced a part of sis¢’ 
and not of rises. (All of the nonads after the sixth were composed entirely 
of supplied zeros both for N and qise.) The last figures of siss’ were copied 
on the third line from the bottom of the fifth strip so that by forming the 
difference between the second and third lines the required value of rise 
came out on the bottom line. This remainder was found on Aug. 11, 1944, 
to be 118 57508 80382 71696 98184 73569 85091 23773 18030 92037. Since 
this residue is not zero it follows that M157 is composite and incidentally 
that My»; still retains the position of being the largest known prime number. 
For every value of k from 8 to 156 the numbers on the three corresponding 
strips were found to satisfy the relation s,’ = N-q, + 7, for each of the 
moduli 10* + 1, 104 + 1 and 10’ + 1. The author desires to announce 
that he has already begun to investigate M¢7. 


1 Archibald, R..C., Scripta Mathematica, 3, 112-119 (1935). 

2 Lehmer, D. H., Jour. London Math. Soc., 9-10, 162-165 (1934-1935). 
3 Powers, R. E., Bull. Amer. Math. Soc., 40, 883 (1934). 

4Lehmer, D. N., Amer. Math. Monthly, 30, 67, 68 (1923). 


ON THE STABILITY OF TWO-DIMENSIONAL PARALLEL FLOWS* 
By C. C. Lin 
GUGGENHEIM LABORATORY, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated August 24, 1944 


1. Introduction.—Heisenberg’s! remarkable contribution to the hydro- 
dynamic stability of two-dimensional parallel flows has not been favorably 
accepted and properly appreciated, because his paper is not completely 
free from obscure points. Nor has the work of Tietjens,? Tollmien® and 
Schlichting‘ been properly estimated. As a result, the theory to account 
for the instability of laminar flow at high Reynolds numbers has become 
very confused, and its further development has been very much retarded. 
Various authors suggest that it is necessary (1) to consider disturbances of 
finite amplitudes, (2) to include the effect of compressibility or even (3) 
to modify the Navier-Stokes equations. The present situation of our 
knowledge may be seen from the general lectures given by G. I. Taylor® 
and J. L. Synge.® 

Recently the present writer carried out some investigations in an at- 
tempt to clarify the situation. The theory of Heisenberg was critically ex- 
amined, somewhat modified and further developed. These developments 
were made with particular emphasis on the general criteria of instability 
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and their underlying physical mechanism, both in an inviscid fluid and in a 
real fluid. According to the general criterion obtained, the plane Poiseuille 
flow and the Blasius flow are both unstable at sufficiently high Reynolds 
numbers.’ Detailed numerical calculations were also worked out in these 
cases. The present article gives a brief account of the essential results ob- 
tained. A detailed treatment will be published elsewhere. 

‘2. The Equation of Orr and Sommerfeld.—In studying the stability of 
two-dimensional parallel flows, we may restrict ourselves to two-dimen- 
sional disturbances, according to an investigation of Squire.6 A two-di- 
mensional periodic disturbance in a field of flow with main velocity w(y) 
parallel to the x-axis may be represented by the stream function y = 
o(y)e**—™ , and the linearized differential equation for ¢(y) is 


Pe 
(w — c)(y” — a®’g) — w"y = ee — 2a%ye” + aty). (2.1) 


All the quantities involved are dimensionless: R is the Reynolds number of 
the main flow, and a, c are constants describing the nature of the disturb- 
ance. We shall regard a as real so that 27/a is the wave-length of the dis- 
turbance. The real part of c gives the phase velocity while its imaginary 
part (multiplied by a) gives the index of damping or magnification. The 
function w(y) is regarded as defined for all complex values of y by analytic 
extension. Physical considerations then lead to certain homogeneous 
boundary conditions to be satisfied by ¢(y) at two real points y; and 2 of 
the complex y-plane, which represent the coérdinates of certain layers in 
the field of flow. There are altogether four boundary conditions. Hence, 
we are led to an eigenvalue problem, and a relation of the type 


F(a, c, aR) = 0 or ¢ = c(a, R) (2.2) 


holds, where F is indeed an entire function of the parameters. Since sta- 
bility, neutral stability or instability of the motion correspond to the real 
part of —iac less than, equal to or greater than zero, the problem is to de- 
termine the imaginary part c,(a, R) of c(a, R) for a set of real values of a 
and R. 

The actual investigation depends upon the solution of (2.1) by the 
method of successive approximations. Two methods are used. First, we 
introduce the variable 


n= (y — y0)/e, (2.3) 


where w(yo) = c, € = (aR)~™*, and w’(yo) is positive for real values of c. 
In this way, we obtain a fundamental system of four solutions y(y) in the 
form 








318 PHYSICS: C. C. LIN Proc. N. A. S. 


ey) = Le exi(n), (2.4) 
where x(n) are explicitly expressible in terms of the Hankel functions 
HY [2/s(ia0n)*"], ao = {w'(yo)}/*. Indeed, the series (2.4) can be 
shown to be convergent for the range of variables and parameters concerned. 
Secondly, we may put ‘ 


(9) = D (aR)"eG), (2.5) 


which gives us two asymptotic solutions. The initial approximation of (2.5) 
satisfies the inviscid equation 


(w — (eo — ate) — w'e = 0, (2.6) 


which can be solved by expanding ¢ in power series of a®, convergent for 
all values of a*. Two other asymptotic solutions are obtained in the form 


o(y) = exp{ fgdy}, — g = (aR)7go + 1 + (aR) 2g2 + .... (2.7) 


Both methods of solution are essentially due to Heisenberg,’ and Toll- 
mien followed the same way. Heisenberg also carried out some approxi- 
mate calculation for the case of channel flow, while Tollmien studied the 
case of the boundary layer. Criticism of Heisenberg’s work is usually as- 
sociated with the multiple-valued nature of the asymptotic solutions and 
the convergence of the power series in a? used for solving (2.6). Detailed 
investigations justify Heisenberg’s treatment. In fact, the asymptotic 
solutions hold only for certain regions of the y-plane to be determined by 
comparison with the asymptotic expansions of (2.4). Indeed, from the 
properties of the Hankel functions, it can be easily shown that under the 
restriction 


_ = < arg (aon) < (2.8) 


T 
6’ 
the asymptotic expansions do not suffer any jump (which may be called 
the Stokes phenomenon or the Uebergangs-substitution). Fortunately, 
such a restriction does not exclude the existence of a connected region en- 
closing the points y, and 2 where our boundary conditions are to be satis- 
fied. The eigenvalue problem can therefore be conveniently treated with 
the help of the asymptotic solutions. However, the region of validity of the 
asymptotic expressions does not always include all the points of the real 
segment (y:, ye). As a matter of fact, the asymptotic solutions hold 
throughout the real interval (1, ye) only when c; > 0. Forc,; < 0, there are 
two points on the real axis where the Stokes jump of the asymptotic solu- 
tions occur. These two points coalesce into one critical point when c; = 0. 
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Physically, the failure of the asymptotic solutions at some points on the 
real axis (and therefore the failure of y as a first approximation) in the 
case of neutral (c; = 0) and damped oscillations (c; <0) indicates the pres- 
ence of the effect of viscosity at the corresponding layers inside the fluid. 
This was noted by several authors in the cases of neutral stability. How- 
ever, the existence of two inner viscous layers for the damped disturbance 
does not seem to have been noticed before. As we shall see in Section 3, it 
has an important consequence upon the well-known criteria of Rayleigh’ 
and Tollmien.° : 

3. General Criteria of Stability and the Curve of Neutral Stability—In the 
limit of infinite Reynolds numbers, it can be shown by careful mathematical 
investigation that the problem can be treated by finding a solution of (2.6) 
satisfying two boundary conditions. Rayleigh’ and Tollmien” have shown 
that the necessary and sufficient condition for the existence of ‘neutral or 
self-excited disturbances is that the velocity curve has a point of inflection 
in the interval (1, ye). The necessary condition is general, while the suf- 
ficient condition has been proved by Tollmien only for velocity distribu- 
tions of the symmetrical type or of the boundary-layer type. In fact, an 
example has been shown by the present writer where a point of inflexion 
exists in the velocity curve, but a neutral disturbance does not exist. 

The authors mentioned above also concluded that damped disturbances 
are excluded if a point of inflexion does not exist. However, their proof as- 
sumes the analytical nature of the solution ¢(y) along the real axis, 
which is true for self-excited disturbances, but not for damped disturbances. 
Previous authors regarded a damped disturbance as the complex conjugate 
of a self-excited disturbance. This is in contradiction with the above-men- 
tioned drastic difference of nature of those solutions along the real axis. 
Careful investigation reveals that the damped or self-excited nature of a 
solution is not changed by taking. the complex conjugate. Therefore 
damped disturbances are not excluded by the condition w"(y) ¥ 0 (at least not 
according to Rayleigh’s proof). 

Tollmien’s proof to establish the existence of a point of inflexion in the 
velocity curve as sufficient to insure instability is rather cumbersome, and 
he has to assume that w’’’ does not vanish at the point of inflexion. These 
points have now been improved by a more careful consideration of the 
analytical nature of (2.6). 

Having settled the question of hydrodynamic stability at infinite Rey- 
nolds numbers, the author studied the stability in a real fluid by trying to 
find out how the foregoing results should be modified for large but finite 
values of the Reynolds number. One such result has already been obtained 
by Heisenberg,' whose conclusion is as follows. If a velocity distribution 
allows an inviscid neutral disturbance with finite wave-length and non- 
vanishing phase velocity, the disturbance with the same wave-length is un- 
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stable in the real fluid when the Reynolds number is sufficiently large. In 
the present investigation, the proof has been formulated in a slightly more 
satisfactory manner. 

Since a point of inflection in the velocity. curve is sufficient to insure in- 
stability only for velocity distributions of the symmetrical and boundary- 
layer types, the author has also made more attempts in these cases. _ In- 
deed, it has been shown that in these cases the flow is always unstable for 
sufficiently large Reynolds numbers, whether the velocity curve has a point 
of inflection or not. The curve of neutral stability 


ca, R) = 0 (3.1) 


may belong to either of the types shown in figure 1. When the velocity 
curve has no point of inflection, the two asymptotic branches of the curve 
have the coxnmon asymptote a = 0 (Fig. 1 a). When there is a point of 
inflection, one branch has the asymptote a = 0, while the other has the 
asymptote a = a, > 0 (Fig. 1 5). These results are in agreement with 
those of Rayleigh, Tollmien and Heisenberg cited above. 


“ a 

















FIGURE 1 (a) FIGURE 1 (b) 
General nature of the curve of neutral stability. The dotted curve is a 


curve of stability given by Synge. 


Simple general formulae have also been derived to express the exact 
forms of those asymptotic branches in terms of the velocity curve. The 
fact that these two asymptotic branches join together to give a maximum 
a and a minimum R (instead of going to infinite a) can be inferred froma 
criterion of Synge.® It states that there is always stability (c, < 0) if 

(gR)? < (20? + 1)(404+ 1)/a?, g= max. |w'(y)|. (3.2) 


msysy2 


Indeed, the minimum value of the Reynolds number on the neutral curve 
marks the beginning of instability and is therefore very important. Ac- 
cordingly, the following approximate formulae have been derived for the 
evaluation of this minimum value: 
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Pio ’ V2 1/2 
R= se 1 = f wdy t , for symmetrical profiles, (3.3) 
v1 / 
R = 21, tor boundary fil 3.4 
7 cary or boundary-layer profiles, ( A) 


where w; is the value of w’(y) at w = 0, and c is given by 


__ 2") 
* {w'(y)}? 


In these dimensionless expressions, the velocities are referred to the maxi- 
mum velocity and the lengths to the distance from w = 0 tow = 1. Com- 
plete numerical calculation of the curve ¢;(a, R) = 0 has been carried out 
for the plane Poiseuille flow and the Blasius flow. The results are in gen- 
eral agreement with figure 1 (a). For the channel, the critical Reynolds 
number is found to be 16,000 as defined in terms of the maximum velocity 
and the width of the channel; for the boundary layer, it is found to be 400 
as defined in terms of the free-stream velocity and the displacement thick- 
ness. The result in the latter case is in better agreement with Tollmien’s. 
original calculation? than with Schlichting’s later one.‘ The maximum 
value of a is even larger than Tollmien’s value, while Schlichting’s value is 
smaller. 

4. Physical Interpretations and Concluding Discussions—We shall now 
try to understand the physical mechanism underlying the above results. 
Generally speaking, the investigation in an inviscid fluid brings out the 
réle of pressure and inertial forces. The subsequent investigation of viscous 
forces will then settle the stability problem of a real fluid. 

Equation (2.1) is the equation of vorticity. This suggests that the con- 
dition w” = 0 in Rayleigh-Tollmien’s criterion means an extremum of vor- 
ticity of the main flow w(y). In fact, the physical mechanism can be 
understood by following this line of thought. Let us regard the field of flow 
as vorticity distribution in parallel layers, and imagine a disturbance con- 
sisting of the interchange of two fluid filaments in two different layers. 
These fluid eléments will retain their original vorticities, and will therefore 
be present in layers with vorticities different from their own. It can be 
shown that these filaments will then be pushed back to their original layers 
if the vorticity gradient of the main flow has the same sign.'? The flow is 
therefore stable. Indeed, if we consider an excess concentrated distribution 
of vorticity ¢(x, y) (a ‘‘vortex’’) upon a small region of a given field of vor- 
ticity, fo(y), then the vortex is accelerated with an average acceleration 


1 
atest - f v(x, y)fo' (y)dxdy (4.1) 


= 0.6 at w(y) = ¢. (3.5) 
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in the direction of the positive y-axis, where I is the total strength of the 
vortex, v(x, y) is the y-component of velocity due to it, and the integration 
is taken over the whole field of flow. A consideration of the sign of I and 
fo’(y) will then substantiate the above statements. Thus, the stability in 
an inviscid fluid as governed by inertial and pressure forces can be discussed 
in terms of the gradient of vorticity of the main flow. 

Another interpretation of the significance of the gradient of vorticity is 
to be found from Kelvin’s ‘‘cat-eye’’ picture of a flow with neutral dis- 
turbance.’* If we consider the vorticity of fluid elements along a closed 
stream line in that picture, we see that such a flow pattern is not possible 
without diffusion of vorticity by viscous forces, if the gradient of vorticity 
is finite at the critical layer where w = c. 

This consideration immediately leads to a simple method of visualizing 
the extent of viscous forces in controlling the stability of fluid motion. If 
we consider the diffusion of vorticity from the critical layer for a period T 
of a neutral disturbance, the effective distance covered’ would be of the 
arder of (vT)” i being the kinematical viscosity. The ratio of this dis- 
tance to that between the critical layer and the solid boundary can be 
easily shown to be 


Qa 
where z = — aon; is directly related to the argument of the Hankel functions 


of Section 2, m being the value of 7 at y = y,. This parameter s has the 
value 0.7 at infinite Reynolds number on the lower branch of the neutral 
curves of figure 1. It decreases to about 0.5 at minimum Reynolds number 
. and further decreases to zero as R becomes infinite along the upper branch. 
Referring to those figures, we note that for small values of s (small viscos- 
ity), the effect of viscosity is essentially destabilizing, since an increase of 
Reynolds number gives more stability (cf. Heisenberg’s criterion of Section 
3). For large values of s (large viscosity), the opposite is true. The physi- 
cal mechanism can be described in the following way. It is known that the 
destabilizing effect is caused by phase shifts of the disturbance, which tends 
to build up its shear —pu’v’ (u’, v’ being the components of velocity of the 
disturbance, and p the density of the fluid). This, in turn, transfers the 
energy of mean flow into that of small oscillations. The stabilizing effect 
is due to dissipation. Now, for small values of s, we may think of the flow 
as having two thin viscous layers, one at the wall, the other at the layer 
where w = c; the effect of dissipation is relatively unimportant, and the 
resultant influence of viscous forces is destabilizing. When s exceeds 0.5, 
these two layers may be regarded as having joined each other into one vis- 
cous region; the effect of dissipation becomes important, and the resultant 
influence is stabilizing. 
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The viscous forces upset the decisive nature of an extremum of vorticity 
for the determination of the stability of the flow of a real fluid, as one can 
see from a comparison of figure 1 (a) (no extremum) and figure 1 (6) (with 
extremum). Thus, there is no drastic change in the stability characteris- 
tics of a boundary layer as the pressure gradient changes from a slightly 
positive value to a slightly negative value (as one would expect if one tries 
to apply Rayleigh-Tollmien’s criterion to a real fluid). However, the ap- 
pearance of the point of inflection in a velocity curve does indicate an in- 
, crease of instability. It can be shown, by using (3.3), (3.4) and (3.5), that 
the critical Reynolds number will be decreased. Also, for a given Rey- 
nolds number, the range of wave-lengths of unstable disturbances will be 
increased. For a boundary layer with adverse pressure gradient, a further 
indication of increasing instability is given by the decrease of w;. Indeed, 
as w; — 0, the approximate formulae (3.4), (3.5) givec—>1,R—0. This 
is, of course, an extrapolation of those formulae, but the conclusion stands 
in qualitative agreement with the usual notion that the flow would be vio- 
lently unstable in such a case (separation). 

5. Transition to Turbulence.—Finally, we want to touch briefly the ques- 
tion of transition to turbulence. A linearized theory of hydrodynamic 
stability is of course not quite adequate to account for such a phenomenon. 
But if we allow our theory to be extrapolated to the beginning of non-linear 
effects, and assume the subsequent development into turbulence to be very 
rapid, there is the possibility of determining approximately the Reynolds 
number of transition. For flow in boundary layers, it is interesting to note 
that the disturbance in the free stream has an isotropic nature in two di- 
mensions. There may be therefore the possibility of introducing the re- 
sults of the theory of instability in boundary layers into Taylor’s theory of 
the transition Reynolds number,’ which appears to be quite successful 
when the transition is caused by the turbulence in the free stream. 


The author wishes to express his sincere gratitude to Professor Theodore 
von Karman for suggesting this problem and for his invaluable help 
throughout this work; to Professor Clark B. Millikan for helpful sugges- 
tions and for using some of his unpublished notes, from which much inspira- 
tion and suggestion has been derived; to Professors P. S. Epstein, H. Bate- 
man, Dr. Hans W. Liepmann and several of his friends for many valuable 
discussions and suggestions. 


* The present article is a brief report of some results contained in the author’s disser- 
tation approved by the California Institute of Technology for the degree of Doctor of 
Philosophy. 
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FLAT SPACE-TIME AND GRAVITATION 
By GreorceE D. BIRKHOFF 
DEPARTMENT OF MATHEMATICS, HARVARD UNIVERSITY 
Communicated September 7, 1944 


If one admits that physical events take place in a 4-dimensional space- 
time continuum (an idea abandoned in current quantum-mechanical the- 
ory) there are three interesting possibilities: classical space and time; 
flat or electromagnetic space-time; curved space-time. The appropriate 
corresponding mathematical languages are, respectively, those of 3-vec- 
tors, 4-vectors and 4-tensors. 

In a certain sense the flat space-time, characteristic of the so-called special 
theory of relativity, is just as absolute as classical space and time, since the 
coérdinates ¢, x, y, require exactly 10 arbitrary constants for their com- 
plete specification in both cases. But, in the framework of flat space-time, 
- the fundamental electromagnetic equations of Maxwell and Lorentz lose 
the artificiality which they possess in classical space and time. 


The initial attempts to incorporate gravitational phenomena in flat 
space-time were not satisfactory. Einstein turned to the curved space- 
time suggested by his principle of equivalence, and so constructed his gen- 
eral theory of relativity. The initial predictions, based on this celebrated 
theory of gravitation, were brilliantly confirmed. However, the theory has 
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not led to any further applications and, because of its complicated mathe- 
matical character, seems to be essentially unworkable. Thus curved space- 
time has come to be regarded by many as an auxiliary construct (Larmor) 
rather than as a physical reality. 

In my opinion, the failure of the early attempts by Nordstrém ,and 
others to develop a theory of gravitation in flat space-time is to be at- 
tributed to the fact that a fundamental theoretic requirement was over- 
looked, namely, that the disturbance velocity in matter must be that of 
light. 

With this requirement in mind, I have recently been led to a very simple 
theory of gravitation in flat space-time, concordant with all known gravi- 
tational phenomena and free of arbitrary constants. This theory was pre- 
sented first in very brief form in 1942 at Tonanzintla, Mexico, and has been 
developed further in a Note in these PRocEEDINGS.! Furthermore, atten- 
tion is to be directed to a Note by A. Barajas,? called forth by a review of 
H. Weyl, and to an article by A. Barajas, C. Graef, M. Sandoval Vallarta 
and myself, taking up the new theory from the physical point of view.* A 
very significant application of the theory to the two-body problem by Graef 
will be published shortly.‘ 

Unfortunately, the foundation of the theory has not so far been ade- 
quately presented in its philosophic, postulational and mathematical as- 
pects. My colleague, Professor Barajas, and I are planning to publish an 
extensive article remedying this deficiency. 

The aim of the present Note is to present briefly these foundational con- 
siderations as I see them. It is especially necessary to do so in order to 
avoid further misunderstanding of the new theory. For example, Weyl 
says very recently,’ referring to my theory: ‘Their [i.e., ‘the field equa- 
tions’] most general static centrally symmetric solution involves 3 arbitrary 
constants a,b,/.... From the present standpoint this is a serious disadvan- 
tage of B.”’ His assertion is wrong since the general exact solution for the 
gravitational potentials h,; is 


m 
hy = by y 


where 7 stands for radial distance, 6, is the familiar Kronecker 4, 
and the mass m is the single arbitrary constant which enters. This 
exact solution plays in my theory a réle analogous to that of the 
Schwarzschild solution in the theory of Einstein. Weyl has overlooked the 
salient fact that the central body is composed of the basic ‘‘perfect fluid.” 

The proposed theory of gravitation in flat space-time may be character- 
ized as follows in its fundamental features: 

1. Inthe 4-dimensional framework of flat space-time, matter is regarded 
as the occupant of certain’ tubular regions made up of the worldlines of 
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identifiable points. Point-particles are abandoned once and for all, except 
as a limiting possibility. 

Thus there is a duality of matter and space-time in my theory, whereas 
the monistic concept of space-time conditioned by matter prevails in Ein- 
stein’s theory. Whichever point of view is destined to figure in ultimate 
physical theory, it seems the part of obvious common sense to explore both 
possibilities fully. 

The simplest available form of matter in flat space-time is that charac- 
terized by a certain stream 4-vector v' of space type, i.e., with 


p? = (0')? — (0°)? — (0)? — (9)? > 0; 


here p is the scalar length of the 4-vector v'. Or, alternatively, we may 
think of matter of this kind as characterized by a density p and a velocity 
4-vector u’, where v’ = pu’ and 


(ut)? — (u?)? — (u%)? — (u4)? = 1. 


The principle of local causation is taken to hold for an isolated portion 
of this kind of matter, as embodied in the differential equations: 


a Pane eee 
Bm Pv. 0 B.S) 


where ? = x1, x = x*, y = x’, g = x‘, and where F' are taken to be rational 
and integral in the partial derivatives involved. These four differential 
equations assert that the time rates of change of the components of the 
stream vector are functions of these components and their space rates of 
change, being rational and integral in the latter. 

The requirements of the underlying language of 4-vectors, based on the 
Lorentz group, then indicate as the only available possibility, when referred 
to instantaneous rest codrdinates (v! = p, v? = v* = vt = 0): 


ov! a eae 
— FS + 22 + BY 
Ov ov! Te 
ar = Ca (i = 2, 3, 4). 


These equations involve a pair of arbitrary functions F(p), G(p). 

By appropriate normalization of the scalar p (that is by multiplicative 
modification of the 4-vector v‘) the preceding equations may be expressed 
in the normal form 


oT oa 


sear Sie ij 
oo 0 (T’ = v'v’ — p(p)g’’) 





where g!! = —g?? = —g33 = —g# = 1, 7 = Ofori+j. Here there ap- 
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pears only one arbitrary function p(p). This normalized density p is de- 
termined up to a unit of magnitude, while the scalar p(p) is determined 
up to an additive constant. 

2. Itis granted that free equilibrium is possible for a certain equilibrium 
density, po >0, i.e., that a possible state is v' = w' provided only that we 
have 


(vo)? ae (v0)? des (v9)? eat (vo4)? = po’. 
Because of the thoroughgoing analogy of the equations obtained with those 
of a homogeneous adiabatic fluid in classical physics, it is natural to assume 
that along a free boundary we have always p = po, and furthermore that at 
collision p takes on the same value on both sides of the common boundary 
until separation occurs for p = po subsequently. 

At this stage the behavior of a collection of freely moving portions of this 
“fluid” has been completely specified, whether or not collisions occur. The 
equations involved present a more familiar aspect if the velocities u* are 
introduced, with 


T# = pu‘u’ — p(p)g”. 


This type of equations has always been taken to be appropriate for a gen- 
eral homogeneous adiabatic fluid in flat space-time. The symmetric tensor 
T* has been called the energy tensor; p, the density; and (p), the pres- 
sure. 

3. Such a fluid has the property that a certain divergence vanishes: 


Se Sie v*) = 0 


This ensures that the 3-dimensional integral 
S pe~ S40/? dy 


over the rest-volume is invariable.® 
Keeping in mind the hydrodynamic analogy, it appears to be absolutely 
essential to suppose that this divergence vanishes under all conditions. 
Otherwise the fluid might undergo a full cyclic return to a set of initial 
velocities without a cyclic return of densities, such as is always observed 
to occur with ordinary matter. This requirement means that we have al- 
ways 
1,2 
6 Ox 
4. There is, however, a fundamental theoretic difficulty in the theoretic 


employment of the general adiabatic fluid. In fact if two portions of the 
fluid collide with oppositely directed velocities nearly equal to that of light 


=f). 
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(1)’ the equations of motion will break down if the disturbance velocity of 
the fluid is less than that of light. Since it is physically inadmissible that 
this velocity v (relative to rest codrdinates) exceed that of light, we are led 
to require that this velocity, namely, 


v = Vdp/(dp — dp), 


equals that of light, so that at all densities dp/dp = 1/2. 

Thus we obtain by integration the only physically allowable constitu- 
tive equation p = p/2. The corresponding special form of the general adi- 
abatic fluid is called the perfect fluid. 

As far as this determination of T” is concerned, we might have written 
more generally p = '/2p + c, and so have obtained an additional term of 
the form —cg” in JT”. This modification would not affect the equations of 
motion, however. Reason will be given later on for the special choice made 
of c = 0 in fixing completely the energy tensor 7”. 

For the perfect fluid the invariable integral reduces to the simple form 


SV pdv 


where dv is the 3-dimensional rest volume. 

The perfect fluid may be looked upon as the counterpart in flat space- 
time of the homogeneous adiabatic incompressible fluid in classical space 
and time, which has infinite disturbance velocity. Physically speaking, the 
perfect fluid is very nearly incompressible and thus possesses very nearly 
invariable mass J pdv. 

If electricity of density o be attached to the perfect fluid, the ratio 
a/ +/p, called the substance coefficient, remains forever constant along any 
worldline. 

In what follows the perfect fluid is regarded as the single primordial form 
of matter. 

5. Suppose now that the perfect fluid, with energy tensor 


T4 = vv’ — 1/,(v,0%)g (1) 
is not free but is subject to body forces. Formally, the suggestion is obvious 
that we define a force vector f’ by means of 
a 
Ox 





r (2) 


where because of the postulated invariance of the integral 
S V pedo (3) 


we have 


f.0° =0 (4) 
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Thus the force vector is required to be identically orthogonal to the 
stream (or velocity) vector. It may be recalled that the acceleration vector 
a‘ along any worldline has the same property. 

In view of the identity (4) it is clear that f‘ cannot be independent of the 
vector v‘. For the case of electrically charged matter f‘ is known to be 
linear homogeneous in the 4-vector v‘ and identically orthogonal to it. For 
the gravitational forces it is the simplest possible hypothesis to assume that 
in a purely gravitational field f‘ is homogeneous and quadratic in v‘ and 
proportional to p. We write therefore in that case 

ft = poigo*v? (5) 
" where the components of the tensor in = ¢}j are functions of t, x, y, 2 de- 
fined throughout space-time. 

Without going into any detail it is to be stressed that this is the most 
natural assumption which can be made about the manner in which f* de- 
pends on the stream vector v‘. This is especially the case in view of the re- 
versibility of gravitational phenomena in time. 

6. At this stage of our genetic account of the theory under consideration, 
in full accordance with the tradition of the past, it will be supposed to begin 
with that non-gravitational forces, such as electrical forces, need not be 
considered in gravitational problems. This hypothesis is justified by the 
more complete form of the theory given in the Note already alluded to.! 

For the case of a rest system the x, y, z components of the gravitational 
forces take the form 


fe =p, = fy = pe, fe = pois (5’) 
with f, = 0 of course. 
In the corresponding Newtonian situation the force components are: 


ae. nc 
ox’ fu = Pa fr = Px 


where g is the Newtonian potential defined by Poisson’s equation 


2 2 2 
oe 5 + oe = — 4np, or 0 in empty space, 
together with the requirement that g is finite (or vanishes) at infinity. 

By formal analogy one is led to require that gj, are linear homogeneous 
in the first derivatives of the gravitational potential defined by means of an 
appropriate Poisson equation. 

A very simple possibility would be to set up a scalar gravitational poten- 
tial h defined by 

Oh Oh Oh Oh 


_ oa dyt ~ Dat = 8nrT, or 0 in empty space, 





« 
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where we have written 
T= Lag l = p/s 


for the contracted energy tensor. However, with this form of f* it is not 
possible to fulfil the vital condition of orthogonality, embodied in (4). 

The alternative, equally simple, hypothesis from the formal point of view 
is to take as the analog of the Poisson equation 

E ~ os — > _ 2 hiy = 8rT;;, or 0 in empty space, (6) 

with the additional requirement that /,; are finite (or vanish) at infinity. 
The gravitational potential h,; thus introduced is a symmetric tensor of the 
second order. No essential limitation is introduced by use of the special 
numerical factor 87 on the right. 

It is now readily seen that the same condition of orthogonality (4) leads 
to the following uniquely determined form for the gravitational force vec- 


tor f*: 
f, (5 re chat f (7) 


Ox? Ox? 


The theory has now been completely formulated for the limited form in which 
only gravitational forces enter. It may be extended soas to include electrical 
and atomic forces in a natural and interesting manner (see reference 1 and 
Section 11 below). In the limited but important form now under consider- 
ation, with use of the usual tensorial subscript and superscript notation, the 
theory is embodied in the pair of equations 


oT 8 iy (2 a Shey 














ox OxF Ox? 
(8) 
af hy _ or or 0 in empty space 
axa Seseroutia 


where 7“ is given by equation (1). 
At this stage it is clear why it was natural to take the undetermined con- 
stant cin T”% to be 0. Otherwise we should have had 


T% = yy! — 1/,(v,0%) g — cg = p(utul — 1/eg"l) — cgtl 


and we should obtain as the Poisson equation for the limiting case p = 0, 
inside of the fluid, 


ap Oy 





DxtOab ~~ Sr eBen 


which seems inappropriate unless c is 0. 
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On the basis of the theory thus obtained it is found that isolated bodies 
have a static, centrally symmetric distribution, and that a collection of 
spherical bodies move with very high approximation according to the 
Newtonian law of gravitation, relative to a convenient rest system. 

7. Let us now turn to the case of a static, centrally symmetric distribu- 
tion of the perfect fluid. It is not difficult to determine the radial distri- 
bution of the density p, but for our present purposes it suffices to observe 
that we have precisely (7, radial distance), 


Ti == 1/5p(r) 54. 
Our extended Poisson equation takes the form 


Ohy , Why ‘ Ohi 
Ox? Oy? oz? 











= — 4rp(r)é,,, or 0 in empty space. 


Here, of course, the boundary of the (spherical) distribution occurs for some 
value 7o of ry where p(7o) = po. Thus we find as the exact solution for the cen- 
trally symmetric case 


m 
hi = — by: 
or oe 


in empty space (r > 79), involving, as the single arbitrary constant, the 
mass m of the fluid. 

In this way we obtain the gravitational potentials around a static, 
centrally symmetric body like the Sun. These potentials are not observ- 
ably affected by random atomic motion. 

8. Now consider a comparatively small approximate sphere of the per- 
fect fluid, forming in the limit a kind of ideal particle of mass 0. First, we 
have essentially p = po throughout, so that we may write 


or SS Se eo 
ont Oye Oe 
with high approximation, since (4) yields with similar approximation 
ou" 
ox 
But we have 
Out dui dx! 
i —— = —— - —— ds? = dt? — dx? — dy* — dz? 
“Ox ds” ds? ( i aE 


along the worldline of the ideal particle, so that we may write 
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oT™ dx «= OT™ = dy OT ds 
ox® «ds? -. * Bis 








ds?’ Ox ds? 


Thus there are obtained the differential equations of motion for a com- 
paratively small body which moves in the field of a larger central body. 
These lead to essentially the same result for the perihelial advance of a 
planet and for the bending of a light photon in the field of the Sun as does 
the theory of Einstein. 

Furthermore, if we assume that the Planck formula 


E=h 


determines the frequency of radiation v, the result of Einstein for the 
“red shift” in light reaching the Earth from the Sun is obtained. Inas- 
much as the precise mechanism of radiation is unknown, it seems more cor- 
rect to employ this basic formula, than to give an explanation in which 
the light-carrying photon plays no réle, such as is afforded by the Einstein 
theory. 

9. A real test of the availability of the new theory in other directions 
is afforded by the problem of two or more bodies. 

As a first approximation to this problem, it is natural to consider the 
limiting case of » “‘ideal particles’ of masses m, m2... Mn, respectively, 
obtained by taking the equilibrium density pp to be very large. It is 
clear that in the neighborhood of each particle P; of mass m, the corre- 
sponding gravitational potentials should have a principal part 

m 
ni 


in instantaneous rest coérdinates. 


Graef has already shown (see references 3 and 4) that the calculations 
involved can be effectively carried out in the case of two bodies of masses 
m, and mz. Presumably his method can be extended to the case of more 
than two bodies. Furthermore it should be possible to investigate in a 
similar spirit all of the fine-structure corrections to the Newtonian theory 
which lie within the limits of observation. 

10. In order to generalize the theory so as to admit cosmological terms 
one has only to write the Poisson equation in the form 

E 22 2 2 


AYE = Ox? — Sy? _ oN = 8rT i; -+- Kgy 


where K is the small cosmological constant. This really means that we 
allow a form of energy tensor 7” containing the term — cg” previously indi- 
cated, with c = — K/8z small but not 0. 

The conditions at © have then clearly to be lightened to the form that 
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hi; becomes regularly infinite to the second order at infinity, in which case 
it is concluded that for a uniquely determined space-time origin we have 


a 
hey = hig + gl = Gh an te By 


where h;; satisfy the previously indicated form of the Poisson equation and 
boundary conditions. In this case there is obtained an expanding uni- 
verse, about the space-time origin ¢ = x = y = z = O in the flat space-time 
under consideration. 

11. In the general case there is an electromagnetic 4-potential ¢, satis- 
fying Maxwell’s equations, and an atomic potential y constant along every 
worldline. This leads to the complete theory with (covariant) force com- 


ponents: 
fe) i fe) a\a@ Oia oh a 
fi= 00 + o( . oe + o( rs oe) PE 





where the terms corresponding to atomic, electric and gravitational forces 
are homogeneous of degrees 0, 1 and 2, respectively, in the velocity vector, 
and are linear homogeneous in the first partial derivatives of the corre- 
sponding potentials y, g‘ and h”’, respectively. 

It will require further mathematical investigation in order to determine 
the serviceability of this conjectural theory in the domain of atomic phys- 
ics. I have previously indicated how an equation much like the Schrédinger 
wave equation may be obtained on the basis of the atomic potential y.® 
However, my attempt was based on a background of curved space-time, 
and I had not then realized that the potentials must all be of zero dimensions, 
so that I used dy/d«* in place of pdy/dx‘. 

12. Professors Barajas, Graef, Sandoval Vallarta and I are examining 
further these and other physical problems in the light of the new theory. 
Meanwhile, it seems clear that the theory promises well and deserves 
careful study because of its striking simplicity, completeness and mathe- 
matical consistency. Furthermore, as will be seen from what precedes, 
the theory is independent of all ideas of curved space-time and of the 
corresponding Einstein theory. 

No doubt, in view of the substantial successes of the Einstein theory, it 
is worth while to attempt to reflect that theory on to flat space-time, and so 
to obtain a degenerate theory, which in a certain sense is only the shadow 
of ashadow. However, the objections made by Barajas? to the form of de- 
generate theory considered by Weyl’ seem to be substantially valid.° 

As far as I can see, the Einstein principle of equivalence (‘that inertia 
and gravitation are one,’ Weyl, loc. cit.) is at bottom a mathematical 
principle signifying only that certain equations in the Einstein theory are 
linear homogeneous in the.density of matter—a fact just as true of the 
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theory of gravitation here proposed as of the general theory of relativity. 
To me there is only the following mathematical fact in the comparison of 
the basic points of view of the Einstein theory and my own: in his theory 
there is no underlying framework of independent variables #, x, y, z valid 
throughout space-time, such as are present in the theory based on flat 
space-time. The real question is whether or not the theory based on such 
special codrdinates is simpler and more useful for the description and pre- 
diction of the physical facts. This is not a question to be decided by a 
priori considerations. What is required is rather a study of the new theory 
and its physical applications. 


1 Birkhoff, G. D., “Matter, Electricity and Gravitation in Flat Space-Time,’ these 
PROCEEDINGS, 29, 231 (1943). My lecture at Tonanzintla is about to appear under the 
title ‘El Concepto de Tiempo y la Gravitacion,’’ Boletin de la Sociedad Matemdtica 
Méxicana, 1, No. 4 (1944). 

2 Barajas, A., “Birkhoff’s Theory of Gravitation and Einstein’s for Weak Fields,” 
these PROCEEDINGS, 30, 54 (1944). 

3 Barajas, A., Birkhoff, G. D., Graef, C., and Sandoval Vallarta, M., ‘“‘On Birkhoff’s 
New Theory of Gravitation,’ Physical Review, 66, 138 (1944). 

4 In vol. 1, No. 5 (1944) of the Boletin de la Sociedad Matematica Méxicana. 

5 Weyl, H., ““Comparison of a Degenerate Form of Einstein’s with Birkhoff’s Theory of 
Gravitation,” these PROCEEDINGS, 30, 205 (1944). 

6 See for instance, Birkhoff, G. D., Relativity and Modern Physics, Chaps. VII and XI, 
Cambridge, 1923, 1927. 

7 The velocity of light c is 1 since the lightsecond is regarded as the unit of distance. 

8 Cf. two notes in these PROCEEDINGS, 13, 160, 165 (1927). 

In his Note Barajas is more than fair to the “degenerate theory’’, which, strictly 
speaking, is no more usable than is the early Nordstrém theory. See M. Wyman, Math. 
Rev. 5, 218 (1944). Since all of the relativistic theories of gravitation take the classical 
Newtonian theory as prototype, the formal resemblance between them is inevitably 
considerable. This fact is stressed, for example, in my article, Newtonian and Other 
Forms of Gravitational Theory, Scientific Monthly, 58, 49 and 136 (1944). It is to be 
looked upon as the source of the formal resemblance between Einstein’s general theory 
of relativity, based on curved space-time, and my own theory, based on flat space-time, 
which Wey] refers to. 














